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ABSTRACT 


A brief description is given of a new refrigerated expansion chamber apparatus based on an Australian C.S.1.R.O, 


design for measuring ice nuclei concentrations. 


The compatibility of results between five instruments of similar 


construction, a comparison of data obtained with a simple version of the mixing chamber method, and homogeneity 


of rapid expansion measurements at sites 8 and 100 miles apart are investigated. 


Except for uncertainties regarding 


the extrapolation of results to natural clouds, all indications are that, with care, the nucleation level in the atmos- 


phere is capable of objective, compatible measurement with standardized equipment. 


However, a series of meas- 


urements verified the existence of significant differences between the rapid expansion and mixing chamber techniques 


Both methods reflected similar trends during appreciable increases or decreases in nucleation activity. 


A surprisingly 


uniform geographical distribution of aerosols responsible for ice crystal nucleation is suggested in some of the results 


1. INTRODUCTION 


The microphysical processes involved in the water-to-ice 
transformation in supercooled clouds has received increas- 
ing attention in recent years. It is generally accepted 
that one of the critical factors is the concentration of 
certain kinds of comparatively rare particles suspended 
in the atmosphere which serve as nucleating agents. 
This particular class of particles is variously referred to 
in the literature as “freezing,” “sublimation,” “deposi- 
tion,” or “ice”? nuclei, depending on the physical mecha- 
nisms invoked. ‘The more general nomenclature of ‘ice 
nuclei’ is used here in view of the current uncertainties 
as to the exact physical nature of the process. 

Since those particles under consideration normally con- 
stitute only an extremely minute fraction of the total 
atmospheric aerosol population, their detection on a 
routine basis at the present time requires the use of 
refrigerated devices and the production of supercooled 
droplets and ice crystals under what must be inferred as 


_—.. 


' These results were obtained in connection with an investigation sponsored jointly by 
the Nat Foundation (under NSF-GL29) and the U.S. Weather Bureau, 


i science 


highly artificial circumstances in the light of the scale of 
events in natural clouds. 


that the concentration of ice crystals produced in a given 


The a priori assumption is 


volume of air is a measure of the number of ice nuclei 
present which are active at temperatures between some 
threshold value and the particular cloud-chamber tem- 
perature at The 
temperature dependency is extremely critical and some- 


which the observations are obtained. 


what variable; a 2° to 5° C. change results in roughly an 
order of magnitude difference in the number of particles 
activated. 

While the basic assumption mentioned above appears 
valid, it is also an over-simplification in the sense that the 
ice nucleation process is known or suspected to involve 
subtle and delicate interactions of such additional param- 
eters as instantaneous supersaturation levels, nuclei 
size, and droplet volume, along with crystallographic and 
Nor is it yet entirely 


the 


surface properties of the nuclei. 


clear whether these nuclei participate directly in 


condensation process, or are activated mainly by coagu- 
lation the exterior 


droplets, or possibly by 


with surfaces of pre-existing cloud 
deposition of moisture directly 
from the vapor state. 
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As a consequence of the complex nature of the problem 
and the present status of instrumentation in the field, it 
has seemed appropriate to adopt an inquiring attitude on 
the question of the measurability of ice nuclei concentra- 
tions in the free atmosphere. The main purpose of this 
paper is to summarize the results of a series of experiments 
designed to partially assess this problem, and incidentally, 
to describe an instrument which has now been in field 
service in multiple numbers over a sufficient period of 
time to be able to report an acceptable level of mechanical 
integrity and durability. 

Although a modest number of ice nuclei measurements 
has been made in recent vears, it has proved exceedingly 
difficult to evaluate or svnthesize the data. Uncertainties 


exist regarding the compatability of instruments and 
measuring techniques, sampling variations in the atmos- 
phere us thes reflect on the representativeness of such 
observations, as well as regarding the knowledgeable 
interpretation of cold box data as precise and meaningful 
The absence of an absolute 
measuring that the 


reliability of such measurements must be, for the time 


cloud physics measurements. 
technique requires judgments of 
being, based largely on evaluations of the internal con- 
sistency of such data. 

The results presented here uppear to contribute con- 
structive and hitherto missing or sketchy information 
with respect to the confidence assignable to such measure- 
ments bevond a laboratory environment, along with indi- 
cations of the compatibility of two methods most widely 
used to date. However, it must be emphasized that a 
large number of important questions on the interpretation 
of such measurements remains to be answered. 

A discussion of the relationships of these and other 
measurements to meteorological events will be 
presented later. rather 
experimental results and s\ noptic clues pointing toward a 


daily 
Some early and unexpected 
potentially significant natural source of ice nuclei, uti- 
lizing the equipment described here, have been reported 
by Brier and Kline [1], Kline [2], and Battan and Riley [3]. 
Most of the measurements discussed in this paper were 
obtained in and around the Washington, D.C. area in 
conjunction with an observational program now in 
progress at a number of widely separated locations during 
the assistance of a number of 
cooperating organizations. Although these data reflect 


the results of an evaluation of two particular techniques, 


selected periods with 


the much broader objective of assessing the reliability of 


such measurements in the free atmosphere strongly 


motivated this study ; 


2. EQUIPMENT 


An external view of one of the identically constructed 
instruments used in this investigation is shown in figure 1. 
The equipment is an adaptation of the basic design de- 
Warner {4}. The essential components en- 
closed in the console consist of: (1) a 10-liter insulated 
copper cvlinder with the wall and bottom external surfaces 


scribed by 
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cloud chamber (Experimental 


Weather 
Australian Com- 


Figure 1. tefrigerated expansior 
Bigg-Warner Type 


modified version of an original design by the 


developed by U.S Bureau as 


monwealth Scientific and Industrial Research Organization 


surrounded by copper tubing functioning as the evaporate: 
in the refrigeration system, (2) a 1/3 h.p. hermetically 
sealed Freon 22 (b.p.—40.5° ¢ 
fan-cooled condenser, and (3) a 1/4 h.p. rotary air pumy 
of 200 ¢.f.m. maximum capacity for pressurization o! 
the chamber. Orifices flush the chamber 
wall are vented to a sensitive differential pressure gage 
(modified air speed indicator) and the air pump. The 
latter is connected to a standard rubber air hose at the 
rear of the console. The dial of the air speed indicator 
is calibrated to indicate a temperature drop on expansion 


‘.) compressor unit witha 


with inside 


equivalent to the theoretical wet adiabatic value. Th 
wall temperature indicator is a conventional dial-type 
bi-metallic thermometer in thermal contact with thé 
exterior of the copper chamber about 1/3 of the distanee 
from the bottom to the top of the chamber. The ai 
seal is maintained between the Plexiglas lid and th 
chamber by adjustable clamps and a Neoprene O-ring 
Decompression resulting in supercooled fog formatio! 
is accomplished by opening a spring-loaded valve in the 
lid after precooling of each air sample. Check and 
relief valves are installed in the air intake system 
retain pressure between the pump and chamber and t 


prevent over-pressurization resulting in damage to tli 


temperature calibrated pressure gage. The expansio! 
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Example of ice erystals produced in a removable tray containing a supercooled sugar solution and placed on the bottom of the 


cloud chamber. This photograph was obtained with a camera mounted at the top of the chamber. 
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ratios used in normal operation are approximately 1.08, 
1.17, 1.26, and 1.37 for observations at —15°C., —20°C., 
25° C —30° C., respectively, depending on the 
wall temperature. 
Special attention was directed to the problem of achiev- 
The system adopted 


., and 


ing stable wall temperature control. 
for primary use relies entirely on a manually adjustable 
automatic expansion valve. This valve serves in practice 
to meter only a sufficient amount of liquid Freon 22 
to the cooling coils to maintain a constant evaporation 
pressure. Manual adjustment is required to compensate 
for the heat load imposed by ambient temperature con- 
ditions. The pressure differential across the expansion 
valve is monitored by reference to a vacuum pressure 
The 
compressor remains in continuous operation, with this 
control thus eliminating the temperature = in- 
stabilitv problems inherent in the conventional thermo- 
statically controlled on-and-off cycling methods. Pro- 
visions for the latter included, however, 
where the equipment must be used in ambient temper- 
ature conditions below which the refrigerant cannot be 
metered in sufficiently small amounts to prevent damage 


gage located on the front panel of the console. 


system 


system are 


to the compressor due to back-pressure. 

Thermocouple surveys on two units operated at the 
normal wall temperature of —10° C. indicated a uniform 
temperature distribution (within 0.3° C.) over the lower 
2/3 of the chamber with an inversion of the order of 
10°C. extending through the top 4 inches, depending 
on the temperature differential between the wall and 
Although introducing some un- 
sampling volume, a temper- 


outside environment. 
certainty on the “effective” 
ature inversion is an essential feature of the expansion 
technique to provide a reservoir of water vapor for re- 
liable supercooled fog production at low expansion ratios. 
Measurements of temperature recovery times of the wall 
and bottom surfaces following purging with air at 23° C 
were all within the range of 90 to 130 seconds, depending 
on the duration of the flushing evele. 

Three particularly attractive features characterize the 
expansion technique. One is the ease with which the 
sugar solution method developed by Bigg [5] may be used 
for objective detection and counting of the ice crystals 
produced in the entire 10-liter sampling volume, as shown 
in figure 2. Second is the rapidity of sampling over a 
spectrum of temperatures through controlled variations in 
the overpressure; and third, a normal operating wall tem- 
perature of —10° C.to —12° C. which considerably reduces 
the risks of contamination from splintering and fragmen- 
tation of frost particles on refrigerated surfaces over those 
encountered under humid conditions at colder wall temper- 
atures. Glycerin is used as a freezing point depressant 
through careful coating of the walls of the copper chamber 
to prevent frost formation. The present version of the 
equipment was deliberately “over-designed” with respect 
to compressor capacity to provide both quick cooling 


capabilities (under 10 minutes from 30° C. to normal 
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operating levels) and optional use at colder wall tempers. 
tures for test with the mixing chamber method. Super. 
cooled cloud temperatures as low as —-30° C. can }y 
obtained either by the expansion process or by selectiy; 
changes in the wall temperature by means of the adjust. 
able refrigeration expansion valve in the case of the mixing 
chamber method.? 


3. PROCEDURES 


The approach adopted for this investigation was die- 
tated by considerations involving the intended use of th, 
equipment. Most of the data were obtained during thy 
course of a twice-daily observational program with th 
instruments used under conditions representative of field 
operation. In principle, greater precision in some of th 
results should be achievable in laboratory environment 
However, possible influences on the results by differen 
operating procedures, sampling variations, and uncertain- 
ties on the nature of the aerosols characterizing those ir 
the free atmosphere were all considered potentially critical 
factors in evaluating the confidence which could be as. 
signed to such The equipment 
located out-of-doors or in a well ventilated enclosure «| 


measurements. Was 
all times to assure purging with representative air samples 
by several rapid up-and-down movements of the tra 
within the chamber. Flushing with the air pump alon 
was found unreliable due to particle scavenging an 
problems associated with non-uniformity of airflow insid 
the chamber. Replication procedures were used through 
out the measuring program. Efforts were made to obtai 
i minimum of five consecutive runs in the case of measure. 
ments confined to a single supercooled-cloud temperatur 
or a sufficient number of observations at varying temper 
tures to define a mean temperature-spectrum-concentr 
tion curve. of —20° C. was 
adopted for all observations confined to a single value 


A “reference temperature” 


The reproducibility of ice nuclei measurements betwee 
instruments (designated by manufacturer’s serial numb 
in the discussion to follow) was examined by simultaneous 
One of tl 
instruments (No. 3) was retained as a “reference” un 
the test 
Efforts were made to cover a representati\ 


operation of two units located side-by-side. 
throughout series, the choice being entirel 
arbitrary. 
range of ice nuclei concentrations to compare both t! 
Each unit was 


frequently operated by different observers with the 1 


mean value and parallelism in the trends. 


crystal counts being made independently. 

An examination of the areal representativeness of expat: 
sion measurements was made in connection with twi 
daily measurements at two fixed sites. One instrume! 
was located on the roof of an office building approximate! 
130 feet above the street level in downtown Washingto! 


D.C. 


A second unit was operated in a suburban res: 


The term ‘mixing chamber method” is usually applied to techniques in which the¥ 


temperature is maintained at a quasi-steady value. The air samples are allowe: gradu 


to reach thermal equilibrium with the refrigerated surfaces. A supercooled clou 


usually produced by introducing moisture into the chamber from an extern uret 
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dentiil section about 8 miles west of the metropolitan 
area. The downtown observations were normally made 
during the morning between 0900 and 1030 Est and during 
the afternoon between 1500 and 1630 est. The normal 
observation times for the second series were 0530 to 0700 
and 1830 to 2000 Est. 
parative series was obtained between downtown Wash- 
ington, D.C. and the U.S. Coast Guard Indian River 
Inlet Lifeboat Station located approximately 100 miles 
These observations were per- 


A subsequent more limited com- 


due east of Washington. 
formed at predetermined schedules during mid-forenoon 
and mid-afternoon periods and included measurements 
over a Spectrum of temperatures in all instances. 

In the case of the comparative runs between the expan- 
sion and the following 


procedure was adopted: A series of five expansion runs 


mixing chamber techniques, 
at —20° C. with different air samples was made with the 
+ ().5° oft 3 
following the expansion series, the wall 
-20° C. 
of measurements with the mixing chamber method. 
One was 


wall temperature maintained at le? C. 
Immediately 
for a similar number 
Two 


iemperature was cooled to 


coud chambers were used for the latter series. 
operated at a temperature of —5° C. for precooling the 
solution, thus permitting nearly instantaneous 


sugar 


response to any ice crystals produced in the second 
chamber after purging and insertion of the cooled tray. 
Counts with the latter method were made at 1- and 2- 
ninute intervals. Two minutes represented about the 
Maximum upper limit practicable with this method due 
to parasitic wall frost problems, unreliable fog duration 
over longer periods, and merging of the growing crystals. 
Each complete set of sequential measurements in this 
lest series spanned periods — nging from 60 to 120 minutes, 
difficulties with 


These experiments were conducted under summer temper- 


depending on contaminating — frost. 
iture and humidity conditions such that a supercooled 
fog formed spontaneously in each fresh air sample drawn 
20° C. for the 


mixing chamber measurements, thus eliminating a need 


into the chamber when operated at 


for the exhalation method or other more elaborate pro- 
cedures for producing supercooled droplets. The sugar 
—10°C 
and sugar by weight plus a wetting agent to reduce surface 
This 


mixture was modified to a slightly more concentrated 


solution used at . consisted of equal parts of water 


tension to provide uniform coverage of the trays. 


strength to delay excessive supercooling resulting in 
spontaneous freezing of the solution and overlapping of the 
ie crystals which would otherwise grow too rapidly at the 
colder temperature used in the mixing chamber method. 
The severity of the frost problem at the operating temper- 
ature of —20° C. was such that it was usually necessary 
(0 reapply glycerin to the wall surfaces after each run with 
the mixing chamber method. Except for those instances 
where the counts were obviously contaminated by wall 
frost, Because of the extraordinary 
tediousness of the mixing chamber method under humid 


ll data were used. 
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environmental conditions, most of these data were ob- 
Tl 
1e 


crystals by the sugar solution technique in both methods 


tained by the same observer. detection of ice 


‘ 


was considered preferable to the “scintillation technique’”’ 
usually used with mixing chamber method, in order to 


reduce sampling uncertainties. 


4. RESULTS 


(a) Compatibility between instruments._A total of 89 
sets of simultaneous runs using 5 instruments in pairs was 
Although a majority of the runs 
20° C., a 


few series included temperature-spectrum measurements. 


performed on 61 days. 
were confined to replicated observations at 


Four representative examples of the latter runs with two 
units are shown in figure 3. A graphical presentation of 
—20° C 
cate the relationship between the readings of the “refer- 
ence” Of these 
instruments, units No. 4 and 5 displaved the widest differ- 


all series obtained at . is shown in figure 4 to indi- 


unit and each of the other units tested. 


ence with a mean ratio with respect to the reference unit 
of 1.5, a statistically significant departure from unity. 
Unit No. 6 displaved values in somewhat closer agreement, 
the average ratios of all the series with this pair being 1.3. 
The comparative series with unit No. 1 gave a ratio of 0.8, 
a statistically significant discrepancy. Correlation co- 
efficients between the paired series ranged from 0.95 to 
0.99. 
sumed te reflect a systematic bias, as is obvious from 


The differences between units were therefore as- 
inspection of figure 4. 

(b) Comparison of rapid expansion and mising chamber 
methods. 
on 74 days. 
intervals are shown as a scatter diagram in figure 5 with 


This test series included 142 sets of observations 


The results corresponding to the 2-minute 


each data point representing the average of the counts 
The 


mathematically fitted regression has significant curvilin- 


obtained in each sequential set of observations. 


earity and indicates that both methods tend to give about 
the same average number of ice crystals only at higher 
the 
extreme concentrations and counting difficulties introduces 


concentration levels. However, rare occurrence of 
some degree of uncertainty in the location of the regression 
curve in this area. The correlation coefficient between the 
rapid expansion counts and the 2-minute mixing chamber 
values was 0.77. Overall, 2-minute mixing chamber counts 
at —20° ¢ 
the ratios varied appreciably over the concentration range 


An examina- 


‘. were 6.8 times the rapid expansion values, but 


as shown by the regression curve in figure 5. 
tion of the trends in mean values between consecutive sets 
of observations revealed that when the rapid expansion 
counts increased or decreased by a factor of 2 (70 cases), 
the trends in the mixing chamber measurements were in 
the same direction in 83 percent of the cases. Changes 
by a factor of 4 or greater (22 cases) resulted in both 
techniques reflecting parallel changes in all instances. In 
un analysis of variance, both techniques showed highly 


significant differences between the means of sets of obser- 
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TEMPERATURE °C 


Figure 3 


chambers a) January 5, 1960; (b 


vations, indicating that the variations of nuclei counts with 
time are real and much larger than the experimental 
“error.” 

If the counts obtained by the rapid expansion method 
are compared with the 1l-minute mixing chamber counts 
(not given here) a correlation of 0.75 results, which is not 
much different than the correlation of 0.77 between the 
expansion chamber counts and the 2-minute interval 
The correlation between the 1-minute interval 
0.84. This 
71 percent) of the 


counts. 


counts and 2-minute interval counts was 


indicates that a large portion (0.84? 
Variation in the counts at 2-minute intervals can be ex- 
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Examples of nuclei temperature-spectrum concentrations obtained during simultaneous operation of two expansiot 
January 9, 1960; (¢ 


January 13, 1960; (d) February 11, 1960 


plained in terms of values at 1 minute. When the counts 
were relatively high at the end of the 2-minute interval 
there was a slight but statistically significant tenden 
(r=0.23) for a higher percentage of the crystals to appea 
during the first minute. 

Another series of measurements data with a differe! 
instrument made over a spectrum of temperatures during 
low concentrations on six different days is shown in figul 
6. These data were obtained by a continuous fog gener 
tion method using a 15-watt light bulb surrounded by 
moist fabric and suspended in the center of the chambet 


a method described by Soulage [6]. The sugar so!utio 
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UNIT NO. 4 (14 SERIES ON /3 DAYS) 
UNIT NO / (26 SERIES ON 1/6 DAYS) 
UNIT NO 6 (22 SERIES ON 16 DAYS) 
UNIT NO 5 (27 SERIES ON /6 DAYS) 


UNITS 


T 


wiTH TES 


QO 


NUMBER 0 


Figure 5.—Results of sequential runs using the rapid expansion and 


te s : , s e ansi ‘assure » st . ‘ ‘ 

Results of simultaneous expansion measurements at mixing chamber techniques at —20° ¢ Each data point is the 
obtained with five units, one of which was used through- average of at least five consecutive samples with each method 
tests as a “control” or “reference” instrument. 


RAPID EXPANSION METHOD 
-——-— MIXING CHAMBER METHOD 





TEMPERATURE °C 


Comparison of rapid expansion and mixing chamber Figure 7.—-Comparision of ice nuclei measurements at two sites 


measurements over a spectrum of temperatures on six different separated by approximately 8 miles 
days. Data points applicable to the same series are indicated by 


simi r geometric configurations 
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was diluted with increasing amounts of glycerin to fore- 
stall spontaneous freezing at colder temperatures. Counts 
were delayed until no additional crystals were observed 
to form. The at 20° (. are considered in 


reasonable agreement with those in figure 5 (mean ratio 


results 
of 5.7 compared with 6.5 for 2-minute values in the same 
indicate that the ratio also 
fog temperature. An 


concentration ranges), but 


varies with the ‘reference’ or 
adjustment of the rapid expansion mean temperature- 
spectrum curve toward 4° to 5° C. warmer values would 
appear to bring the measurements into fair agreement. 
in the concentration range surveyed these results com- 
pare favorably with the conclusions reached from limited 
series of tests with several kinds of instruments as sum- 
marized by Mac ‘ready, Smith, Todd, and Beesmer {7}. 
(c) Areal representativeness of ice nucle: measurements. 
The more extensive series of observations performed to 
evaluate the areal representativeness of ice nuclei measure- 
ments was conducted during January and early February 
1960. Although absolute 
values occurred on a few days, the overall trends were 
which the data are 


considerable differences in 
similar, as indicated in figure 7, in 
presented as 3-day moving averages centered about the 
middle day. There was an apparent tendency for the 
highest counts to occur in the metropolitan area, the 
mean monthly values differing by a factor of 1.65. How- 
ever, it can be shown that this is not a real but an in- 
strumental effect. The units used for this series (Nos. 1 
and 4) 
were compared with unit No. 3, 
When differences in the characteristics of the instruments 
are allowed for, the mean monthly values at the two sites 


have not been operated side-by-side, but both 
as shown in figure 4. 


separated by S miles differ by a factor of 1.26, with the 
lowest counts in the metropolitan area. The correlation 
coefficients in daily mean values and the 3-day moving 
averages shown in figure 7 were 0.71 and 0.85 respectively. 

A further examination of the spatial compatibility 
of ice nuclei measurements was provided by a limited 
series of observations at the U.S. Coast Guard Indian 
River Inlet Lifeboat Station near the Atlantic shoreline 
100 miles east of Washington, D.C. during March 1960. 
These data were obtained primarily for another purpose 
to be reported later. However, this series was coordinated 
with those obtained in downtown Washington, D.C., and 
included 1S simultaneous sets of data at both sites over 
the entire temperature spectrum during similar flow 
These data are summarized in 


regimes on 13) days. 


table | 


TABLE | Comparison of mean values of simultaneous ice nuclei 
observations during similar flow regimes 
No. of ice crystals per 10 liters 
Threshold 
temp Cc 
oo? C °C FC 
Shoreline 17.7 6 76 int) 
Washington, 1.« 17.3 y 119 3, 000 


Estimated 
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The two units used in this series were Nos. 3 and 
4. The latter was used at the Washington, D.C. site. 
and, as shown in figures 3 and 4, tended to indicate the 
higher values. When the ‘correction factor’ of 1.5 js 
applied to the shoreline data it is apparent that the 
counts were similar in the warmer ranges of the temper- 
ature-spectrum concentrations. Unless these results are 
fortuitous, which seems unlikely from close inspection 
of the data, this indicates a surprising uniformity in the 
dispersion of the responsible particles. At the colder 
end of the spectrum, however, the discrepancy of a factor 
30° C. is believed considerably larger 
than This 
that observations in the colder ranges may be particularly 


of more than 6 at 
the probable experimental error. 


suggests 


influenced by pollution from urbanized regions. 


5. DISCUSSION OF RESULTS 


In assessing the outcome of the experiments reported 
here, it is perhaps germane to point out that in some 
respects these tests were rather severe. The measure- 
ments were under environmental conditions representativi 
of field usage. By design, several observers participated 
in the effort and many of the observations were made 
during meteorological conditions in which the reliable us: 
of refrigerated devices is awkward and tedious to a 
extreme. Moreover, the fact that the tests were con- 
ducted in a densely populated area with the ever-presen 
from domestic and other loca 


risks of contamination 


sources conceivably introduced greater Variations that 
might exist elsewhere. 

Judgements on the quality of agreement bet weer 
instruments and techniques depend, of course, on th 
objectives of any measuring program. Day-to-day fluetu- 
ations in ice nuclei concentrations May span up to severs 
greater no 


Clos 


inspection of the rapid expansion data indicates that 


orders of magnitude with factors of 5 or 
uncommon within individual sampling periods. 


fluctuations in time of the order of a factor of 2 or mor 
in the mean values of replicated series were detected b 
all the instruments. The absolute values appear capabl 


of measurement to within at least 25 percent on tl 


average. Differences between instruments (and observ: 
ers) of the order of 50 percent with respect to the ‘refer- 
ence’ unit and up to a faetor of about 2 overall ar 
indicated, however, without special attention to precis 
calibration of the equipment. 

The results of the comparative runs between the rapi 
expansion and mixing chamber methods raise som 
potentially crucial questions on the measurability problem 
Some, perhaps most, of the scatter in the results is attrib- 
utable to uncontrolled factors, particularly with the latt 
technique. These include not only fluctuations in th 
aerosol population over the sampling period, but als 
non-uniformity in supercooled cloud density in_ th 
chamber from day to day, and small differences in t! 
strength of the sugar solution influencing the growth rat 


of the ice crystals, along with the previously mentione 
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probiems with wall frost. Moreover, although the mixing 
chaniber technique is sometimes referred to as a ‘‘constant 
temperature” method, this is far from the case in practice 
due ‘o the instantaneous fog formation with the flushing 
of the chamber and the time required for the temperature 
of fresh air samples to reach that of the wall. This 
extends over an appreciable part of the 2-minute period 
in which ice crystal formation was kept under surveillance 


and varies depending on the temperature differential 
between ambient and operating levels. The mean dis- 
crepancies between the two methods, however, appear 
genuine. A satisfying physical interpretation of these 
differences is exceedingly difficult at this stage because of 
the elusive and contradictory nature of the evidence 
concerning the basic physical-chemical processes involved 
in nucleation phenomena. 

Aside from subtle properties of the nuclei themselves, 
such factors as supercooled droplet size and the time 
dependency of the nucleation process are probable influ- 
The recent work reported by Hoffer |S] points to 
significant influences by soluble substances in depressing 


ences. 
the freezing temperature of droplets. In general, each 
of these factors would operate in the direction of reducing 
the number of ice crystals observed in the rapid expansion 
method because of the smaller droplet sizes and extremely 
short time-intervals at which minimum temperatures are 
maintained compared with those characteristic of mixing 
chamber or similar techniques. 

We are inclined to interpret our evidence as providing 
that either of the 
examined may be used for studying relative variations in 


empirical confirmation techniques 


concentration levels. The rapid expansion method was 
much simpler to use in practice under the prevailing 
experimental conditions. However, some rather elabo- 
rate variant of the mixing chamber method such as 
described by Bigg and Meade [9] may be necessary to 


place such measurements on a more physically realistic 
basis. The alternative would be to calibrate the expan- 
Our 


evidence on the suitability of this approach is not entirely 


sion method against an agreed-upon = standard. 


clear. The data in figure 6 offer considerable promise in 
this regard. However, the convergence tendency in the 
higher concentration ranges noted in connection with the 
results shown in figure 5 hints that a lack of stability in 
any such calibration may be expected. The reason for 
this is not obvious, but the possibility exists that day-to- 
day differences with respect to the nature of the nucleating 
particles and their response to the cloud chamber charac- 
teristics may be critical factors. A similar trend has 
been reported by Soulage [10]. 

There is, of course, the basic question at this point of 
whether the highly artificial conditions (especially super- 
saturation levels) in any refrigerated devices of reasonable 
dimensions provide measurements directly relevant to 
Practically no reliable information 
is available on this point because of the inadequacy of 


natural cloud events. 


Xistiig instrumentation for measuring either nuclei or ice 
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crystal concentrations in natural clouds. Moreover, quali- 
tative observations indicate that natural cloud nucleation 
may often be facilitated by the sedimentation of pre-exist- 
ing ice crystals from higher cloud formations into the su- 
percooled portions of a cloud system and possibly also by 
the splintering and fragmentation of freezing particles. 
Natural nucleation under these circumstances would be 
more efficient than inferred from ‘cold box” 
ments. 


meuasure- 


There remains the difficult problem of choosing a stand- 
ard reference temperature in any compatible measuring 
system, since the mixing chamber and similar methods are 
not readily adaptable to both frequent sampling and tem- 
We arbitrarily adopted 
—20° C. when runs were confined to a single temperature. 


perature spectrum measurements. 


However, as shown in several examples by Kline [2] the 
mean temperature spectrum curves often undergo wide 
excursions in slope and shape so that measurements at any 
given temperature do not necessarily reflect conditions at 

20° C. do 
—15° C. and 


There is, however, a statistical tendency 


other values. That is, high concentrations at 
not always indicate abnormal values at say 
warmer values. 
in this direction, and —20° C. appears to be an adequate 
compromise for some purposes. 


6. CONCLUDING REMARKS 

These results have been presented as a contribution to- 
ward placing some of the measuring problems in better, 
though still incomplete, perspective. Our evidence is that 
there are substantial grounds for believing that compatible 
and representative Measurements can be achieved with 
carefully designed instruments and with attention to pro- 
cedural quality control. At the same time, these results 
emphasize the existence of important instrumental and 
sampling problems requiring consideration in the design of 
future equipment and standardization of measuring tech- 
niques. That the natural aerosols responsible for ice crystal 
nucleation are other than of purely local origin and may be 
widely and fairly uniformly distributed at any given time, 
at least in the region and under the synoptic conditions sur 
veyed, is strongly hinted. This would appear to have in- 
teresting implications in the cloud physics field meriting 
further investigation. 
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ABSTRACT 


Low-level, constant volume balloon (tetroon) flights made from Yucea Flat in the Nevada Proving Grounds of the 


Atomic Energy Commission are utilized to yield an estimate of the nature and magnitude of vertical air motions in 


desert areas. 


It is found that during the day the tetroons oscillate as much as 10,000 feet in the vertical with vertical 


velocities often exceeding 5 kt., whereas during the night the vertical oscillations are on the order of 100 feet and ver- 


tical velocities seldom exceed 0.5 kt. 


There are indications that during the day, at this particular site, helical cir- 


culation patterns exist with axes parallel to the north-south oriented valley floor and with alternating circulation 


sense across the valley, while during the night the small vertical oscillations are more nearly tied to the local topog- 


raphy. 
vertical oscillations associated with solar heating. 


1. INTRODUCTION 


The constant volume balloon (tetroon) flights previously 
reported on by the writers were launched from coastal 
radar sites at Cape Hatteras [1] and Wallops Island [2] 
and most of the tetroon trajectories lay over water. As 
the sea surface presents a fairly homogeneous surface with 
respect to radiative and frictional influences, it is not sur- 
prising that most of these tetroon flights (at least those 
that were reliably positioned) appeared to remain very 
nearly at constant level, the exception being two flights 


which were released within, or close to, a cold frontal 


zone |2|. Because one of the chief uncertainties with 
respect to these superpressured tetroons is the degree to 
which they respond to vertical air motions, it was thought 
desirable to make some tetroon flights in desert regions 
where the vertical air motions would be at a maximum. 
Consequently, during September 1960, 15 tetroon flights 
were made from Yueca Flat in the Nevada Proving 
The tet- 
roons were tracked and positioned by means of the newly 
acquired and M-33 radar of the Weather 
Bureau’s Las Vegas Research Station. This paper pre- 
sents some results obtained from this series of flights, with 


Grounds of the Atomic Energy Commission. 


modified 


emphasis on the nature and magnitude of the vertical 
ur motions indicated by the tetroons. 


2. PREVIOUS WORK WITH NEUTRAL BALLOONS 
Several meteorological groups have used neutral or no- 
lift pibal balloons to estimate vertical air motions near the 


— 


Rese undertaken as portion of programs sponsored by Reactor Development Di- 
¢ Energy Commission, and Air Pollution Division, Public Health Service 


Pure mountain influences on the tetroon trajectories could not be clearly delineated owing to the tremendous 


earth’s surface. Probably the most extensive use of this 
technique was at Oak Ridge, Tenn. between January 1949 
and February 1950, as reported by Gifford [5]. Gifford 
also cites German work along this line, in particular the 
work of Lange [8] in connection with gliding activities. 

Superpressured, constant volume, Mylar balloons (tet- 
roons) offer some advantages with respect to no-lift pibal 
balloons, particularly for the determination of vertical 
air motions. These are: 

1. The tetroon, being for all practical purposes non- 
extensible, will not undergo a change in buoyancy force 
due to variation in balloon volume associated with heating 
effects, as will the pibal balloon. In particular, tetroons 
superpressured to 50-100 mb. have the capability of going 
through sunrise and sunset flight 
altitude, although this capability has not yet been utilized. 


without change in 

2. The tetroon can, within certain limits, be set to 
float at a given altitude, whereas the no-lift pibals must, 
by their nature, be set to float close to the earth’s surface. 

3. The radar-reflective tetroon can be tracked to greater 
distances than can no-lift pibals by double theodolite 
methods, and furthermore tetroon tracking can be ac- 
complished in clouds, at night, ete. 

These advantages make it worthwhile again to carry 
out experiments whose aim is the delineation of 3- 
dimensional air trajectories. In the following, comparison 
is often made between the no-lift pibal data of Gifford 
and the tetroon data obtained at Yucca Flat. 


3. FLIGHT PROCEDURES AND DATA OBTAINED 


Table 1 shows the time of tetroon release, the flight 
equipment utilized, the flight level, duration of flight, 








TABLE 1 









A verage flight 


Dat Ll, altitude 


iunch site 


Launched 
ST (feet) 


1 019-60 )=Radar 1213 6200 
19-60 Radar 1346 4100 
$ %19-0 = Radar 1506 3100 
t 19-40 = =Radar 1749 
19-00 «Radar IS55 1000 
fi 919-40 Radar 2008 160 
7 20-00 ~Radar 0635 1200 
s 20-0) - Radar 093s 
4 O20 40 Yucea Wx. Sta 1204 3500 
1 2-40 =Yueca Wx. Sta 1456 4200) 
1 240 «Yucca Wy. Sta IS13 
12 G20 Radar 1925 
13° 20-00 “Radar 2136 1%) 
14°* 9-21-00 Radar 0702 6300 
i 21-0 Yucca Wx. Sta OSS! ORO 


*Immediately after release, momentarily caught in power line and grounded 
* Towed aloft by radiosonde balloon 


Equipment: T=tetroon, C R=corner reflector, M =radar reflective mesh, L=pib 


Flight termi- 
nation was due to grounding of the tetroon, disappear- 


and other information of possible interest 


ance of the tetroon behind the mountain ridges surround- 
ing Yucca Flat, and (most often) loss of radar contact 
with the tetroon as it approached a mountainside due 
to the appearance of terrain echoes in the radar gate. 
It is seen from table 1 that reasonably long-range tra- 
jectories could be obtained only for flights at relatively 
high level or for flights released upwind from the radar 
site. 

The positioning accuracy of the M-33 radar is not 
known with certainty. It has been estimated [3] that 
the average error in azimuth and elevation angle is 
about 2 mils (0.11 degrees), leading to height and lateral 
displacement errors of about 100 feet 
50,000 feet. The FPS 
at Wallops Island gave errors of less than 10 feet at such 


at 
16 radar utilized 


ranges near 


very accurate 
ranges. However, the loss of accuracy in the M-—33 is 
partially offset by its mobility, simplicity of operation, 
and by the fact that a much smaller crew (2 
handle it. 

As a consequence of the above positioning errors, at 


men) can 


Yuecea Flat the tetroon velocity was determined as an 


uverage over l-minute intervals rather than the 30- 
second intervals utilized on the Wallops Island flights. 
The horizontal velocity data were derived from dial 


readings giving tetroon slant range in vards and azimuth 
angle in mils. However, the digital elevation angle data 
could not be obtained conveniently owing to the position 
of the dial and therefore the tetroon height was scaled 
directly from the plotboard. From the change in tetroon 
slant range, azimuth angle, and height with time, values 
of the longitudinal (V,), vertical 


transverse (VV), and 


(wr) velocity components were obtained. 


4. TETROON TRAJECTORIES 


Figure | shows the horizontal trajectories of the 12 


successful tetroon flights. The tetroons were released 


il lig 
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Tetroon flights from Yucca Flat 






Distance 
tracked 
at altitude 
(n, miles 






Termination due to 


Equipment 





70 10.3 | 422° T, CR Radar lost contact 
is 6.7 | 42” T,2CR Grounded 
58 7.9 | 42° T Y Radar lost contact 
0 00 4 7. 2 = * Radar lost contact 
15 1.9 ie & Re Radar lost contact 
iM 6.1 60” T,2CR, L Radar lost contact 
101 2.8 | 0” T,2CR, M Grounded 
0 0.0 wo T,CR Burst 
sl 15.3 | 0” T,2CR Radar lost contact 
) 13.2 | oO” T,2CR Behind hill 
0 0.0 6” T,2CR,M,I No radar contact 
0 00 @’ T,2CR,M,I1 Caught in power line 
44 2.8 | @” T, CR, M, I Radar lost contact 
52 25.3 | 0” T, CR, M Behind hill 
“4 18.9 6° T,M Behind hill 


from the radar site and from Yucca Weather Station about 
7 miles to the south. The trajectories are based upon 
tetroon positions every minute but for simplicity in the 
figure the positions are indicated only at 5-minute 
intervals. For the sake of clarity, the surface topography 
has also been traced only at 1,000-foot intervals. 

The the 
existence of preferred corridors of escape for air flowing 
northward up Yueea Flat during the midday hours, 
with one corridor extending northward from Whiterock 


few tetroon trajectories available suggest 


Spring (as illustrated by the trajectories of flights 1 and 9 
and the other corridor extending northward from Oak 
the trajectories of flights 14 
A priori, one would not expect that the air 


Spring (as illustrated by 
and 15). 
flow at such heights (3,000-6,000 feet above the earth’s 
surface) would be influenced by terrain during the hours 
of daytime instability, and indeed smoke experiments 
have shown no evidence for such a channeling effect 
The evidence for such an effect from the few tetroon data 
is sufficiently strong, however, to warrant a test. series 
of serial tetroon releases in order to confirm whether or 
not tetroons released from the given area will ever 
traverse the highlands between Whiterock Spring and 
Note that flight 2 


this region, suggesting both an unwillingness for air 


Oak Spring. grounded on a spur in 
parcels to ascend over these highlands and an uncertainty 
with regard to which of the above suggested corridors to 
take. 

Possibly associated with the existence of preferred cor- 
ridors of escape is an intriguing tendency for the azimuth 
of the end point of the tetroon trajectories to shift regular 
in a clockwise sense during the afternoon hours. This 
tendency is shown in figure 2 by the dots. A regressio! 
line, fitted to the afternoon azimuths, indicates an azimuth 
turning of 9 degrees per hour, precisely the pendulum-day 
period at this latitude. It is not certain whether this re 
During times of sea breeze 
turn 


sult is coincidental or not. 


low-level winds obtained from anemometers often 
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Trajectories of the 12 successful tetroon flights at Yueca Flat, Nevada Proving Grounds. 
Dotted segments of trajectories estimated from visual sightings. 
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FIGURE 2 Azimuths of end points of tetroon trajectories (dots) as 
A regression line (solid line 


4 


function of mid-time of flight 
fitted to the afternoon azimuths 


with approximately a pendulum-day period [4] and also 
the Eulerian-type winds obtained from various sites in 
Yueea and Jackass Flats indicate a somewhat similar 
change with time [12]. The tetroon data suggest that 
these wind-direction changes with time extend to quite 
high elevations and may not be purely local in nature. 
Two of the tetroon flights (6 and 13) went southward 
down Yucea Flat, embedded in the nighttime drainage 
flow. Flights 4 and 5 were released early in the evening 
in an attempt to catch the transition from daytime up- 
slope to nighttime drainage flow, but unfortunately both 
flights were lost by the radar as they approached the 
mountainside. In any event the shift might never have 
been detected from these flights bee2use of the probability 
that the height of the tetroons above the earth’s surface 
exceeded the depth of the drainage flow during such an 
early stage of development. In a sense flight 7 detected 
the transition from drainage to upslope flow with the 
flight starting out to the southeast and then curving 
counterclockwise and passing out of the Flat to the north- 
east. In the following the flights are grouped and dis- 
cussed according to whether they took place near midday, 
during early morning or evening, or during the transition 


period. 
TETROON HEIGHT FLUCTUATIONS NEAR MIDDAY 


Figure 3 shows tetroon height as a function of range for 
the five flights (1, 2, 3, 9, 10) released during the early 
afternoon. The essentially north-south height profiles 
have been positioned relative to the radar site in order 
that any tendency for similarity among the height traces 
might be noted. In this case there is only the weakest 
sort of evidence for ascending motion upstream from the 
radar, descending motion immediately downstream from 
the radar, and ascending motion about 30,000 ft. down- 
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stream from the radar. Note that in figure 3, and in 
subsequent figures of this type, the horizontal and vertical 
scales are identical so that the trajectory slope in the 
longitudinal-vertical plane is obtained immediately. This 
slope exceeds 30 degrees for long distances along many of 
This is the order of stream- 
line slope found by glider pilots in mountain-wave regimes 
[13]. 


The magnitude of 1-minute-average vertical velocities 


these afternoon trajectories. 


derived from these five afternoon flights is shown in 
figure 4. Nearly one-fifth of the time the absolute magni- 
tude exceeds 5 kt., with some evidence that in the mean 
the ascending motions are of greater magnitude than the 
descending motions. On several occasions the vertical 
velocity of the tetroon exceeded its horizontal speed. 
Table 2 gives an estimate of the average trajectory 
longitudinal-vertical 


wavelength and = period (in the 


plane) for each afternoon flight. It is emphasized that 
these represent Lagrangian scales and do not necessarily 
vield the length and period appropriate to a streamline 
pattern. While the wavelengths and periods are difficult 
to delineate and some subjectivity is undoubtedly in- 
volved, basically the wavelengths are on the order of 
5-10 mi. and the periods of vertical oscillation are on the 
order of 30-45 min. By comparison, the Oak Ridge 
no-lift pibal data of Gifford suggested wavelengths and 
periods on the order of 5,000 ft. and 15 min. under con- 
ditions of strong instability and light wind speed, and 
wavelengths and periods on the order of 15,000 ft. and 
25 min. under conditions of moderate instability and 
strong wind speed. Thus the period of oscillation in 
the longitudinal-vertical plane is about twice as great 
at Yucca Flat as at Oak Ridge while the wavelengt! 
is about four times as great. 

There exists a ‘‘natural’’ period of air parcel oscillation 
in the vertical [7, p. 132], which period is infinite for a 
dry adiabatic lapse rate and is about 10 minutes for the 
lapse rate (0.7° C./100 meters) usually prevailing in the 
atmosphere. Both the Wallops Island tetroon flights 
and the Brookhaven no-lift pibal flights of Cifford [6 
indicate a period of oscillation in the vertical somewhat 
smaller than 10 minutes for the given stability, in rougl 
agreement with the ‘natural’ period. The 40-minut 
periodicity derived fron the Yucea Flat flights v ould 
correspond to the ‘natural’ period only if the lapse 





TABLE 2.--Lagrangian wavelength and period of vertical oscillation 
for early afternoon tetroon flights from Yucca Flat 


Flight Wavelength Period 
feet minutes 


1 42, 000 

» 26, 000 4 
3 35, 000 # 
4 43, 000 3 
10 51, 000 4 


Average 39, 000 
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FIGURE 3. 


Tetroon height versus distance along trajectory for early-afternoon flights (1, 2, 3, 9, 10 


at Yueea Flat. Dashed lines represent 


interpolated heights. 


rate at Yuecea Flat during the time of the flights was 
within 0.02° ©./100 meters of the dry adiabatic lapse 
Both the Las Vegas and Jackass Flat soundings 
do indeed show that the lapse rate is very close to the 


rate 


dry adiabatic at the time of these flights, with a super- 
adiabatic lapse rate close to the ground followed by an 


approximately dry adiabatic lapse rate up to at least 
10,000 feet. 
fluctuations experienced by 


Thus there is a suspicion that the height 
tetroons and no-lift  pibals 
may be attuned to a period specified by the existing 


lapse rate. 
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1.—Frequency distribution of 1-minute-average vertical 
ties derived from five early-afternoon tetroon flights 


With regard to the fidelity with which tetroons follow 
the vertical air motion, it is of interest to compare the 
height which air parcels might be expected to attain 
during the heat of the afternoon (mixing depth) with the 
height attained by the tetroons. Figure 5 shows early 
morning (0300 pst) radiosonde soundings for Jackass Flat, 
the of Yueca Flat. The 
dashed lines in this figure represent dry adiabats drawn 
through surface temperature readings at Yucea Flat dur- 
and 3 on September 19, flights 


some 20 miles to southwest 


ing the time of flights 1, 2, 
9 and 10 on September 20, and flight 15 on September 21. 





an 
0300 PST Sept 2! 





0300 PST Sept '9 0300 PST Sept 20 











Figure 5.—-Early-morning radiosonde soundings at Jackass Flat, 
Nevada Proving Grounds 


through surface temperature readings at Yueca Flat during times 


Dashed lines represent dry adiabats 


of tetroon flight. 
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TAPLe 3 
mixing de pth for daytime tetroon flights from Yueca Flat 


Comparison of marimum tetroon height and air parcel 


Maximum 
Mixing depth 
tetroon 
feet 
height (feet 


10, loo 


The flat temperature maximum at this locality enables a 
single temperature to represent conditions during more 
than one flight. 
free vertical mixing during the afternoon up to the height 


Presumably, one would anticipate rather 
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where the dry adiabats intersect the sounding cur es, 
Actually, since a superadiabatic lapse rate prevails near 
the earth’s surface in this region during the daytime hours, 
the mixing depth would be somewhat less than that indica- 
ted by the given intersection. Table 3 gives a comparison 
between the mixing depth so obtained and the maxi- 
While there js 


good agreement on some of the flights and poorer agree- 


mum height attained by the tetroons. 


ment on other flights, there is a suggestion that basically 
the tetroons were oscillating in the vertical through uw 
laver of air which one would expect to be in a mixed condi- 
tion. This represents evidence that the tetroons were 
following the vertical air motions to a considerable degree. 
Apparently, the restoring force acting to displace the 
tetroons back to their equilibrium density surface is small 
compared to the drag force acting on the tetroons through 
the agency of large vertical motions. 
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FIGuRE 6 
lines represent interpolated heights. 


Tetroon height versus distance along trajectory for early morning and evening flights (4, 5, 6, 7, 13 


at Yucca Flat. Dashed 


Dotted line indicates surface topography along flight 7. 
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FIGURE 7. 


Tetroon height versus distance along trajectory for morning flights (14 and 15) at Yucea Flat. 


Dashed lines represent inter 


polated heights, dotted lines indicate surface topography along flights, and solid and dashed arrows indicate positions of maximul 


and minimum tetroon height. 
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6. TETROON HEIGHT FLUCTUATIONS AT NIGHT 


Figure 6 shows tetroon height as a function of range for 
the five flights (4, 5, 6, 7, 13) released during the evening 
and early morning hours. Presumably, these flights are 
reasonably representative of flight conditions during night- 
time. The difference in magnitude of tetroon oscillations 
in the vertical at night and near midday is striking 
indeed, with vertical oscillations in the former case on the 
order of 100 ft. and in the latter case on the order of 
1,000-10,000 ft. Furthermore, in contrast to the results 
presented in figure 4, 84 percent of the 1-minute-average 
vertical velocities derived from these five flights were 
within the range —0.5 to 0.5 kt. 

The two flights released latest in the evening (flights 6 
and 13) flew no higher than 300 and 600 ft. above the 
ground, respectively, apparently within an inversion 
surface. It is assumed, therefore, that the vertical oscilla- 
tions experienced by these flights to some extent mirror 
the undulations in the Although 
flight 6 closely followed the contour of minimum height 
in the valley floor whereas flight 13 followed the contour 


inversion surface. 


of 4,080 ft. somewhat to the west, along both flights, in 
general, the tetroon height slowly and regularly increased 
with time suggesting either that the inversion surface 
rises aS One Moves southward down the Flat or that at 
this time in the evening there is an increase in inversion 
height with time. The possibility that the tetroons did 
not reach their equilibrium altitude is made unlikely by 
previous experience concerning the length of time it takes 
tetroons to attain that altitude. The use of tetroons as 
an inversion-surface tracer (in the Los Angeles Basin, for 
example) is immediately suggested. 

Other than the slow increase in tetroon height with 
time noted along flights 6 and 13, little in the way of 
regular height fluctuations was noted on these nighttime 
flights. The rather large, but intermittent height fluctua- 
tions along flight 7 will be considered in relation to moun- 
tain effects in the following section. 


INFLUENCE OF TOPOGRAPHY UPON TETROON 
HEIGHT 


One of the features it was hoped could fairly easily be 
investigated by means of tetroon flights in the Yucea Flat 
How- 


ever, since many of the tetroon flights were confined to 


region was the air flow over mountainous terrain. 


Yucca Flat itself and, moreover, since most of the flights 
were made during the daylight hours when the vertical 
tetroon 


oscillations associated 


outweighed 


with solar heating far 
with 
The only 
flights suitable for consideration here were those released 
early enough in the morning to avoid the bulk of the 
thermal and which 


any vertical oscillations associated 


topography, the results were disappointing. 


influences vet crossed mountain 
ridges of appreciable size (flights 7, 14, and 15). 

The topography along flight 7 
igure 6. This flight 
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has been plotted in 


traversed a mountain pass and 
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Figure 8.—Portions of trajectories of flights 14 and 15 (solid lines 
in the longitudinal plane and their relation to a schematic moun- 
tain ridge. The 
theoretically by Queney for medium-sized mountains. 
line is locus of maximum tetroon heights. 


deduced 
Dotted 


dashed lines represent streamlines 


subsequently disappeared behind a hill for over 10 minutes 
The inter- 
esting thing is that while this flight was stabilized at a 
level sufficiently high to avoid an upcoming mountain 
ridge, upon passing a small hill of no more than 100-foot 
height, the tetroon descended 400 feet to a level below 
the ridge. 


before again being picked up by the radar. 


It is suggested that this tetroon descent was 
induced by descending air motion associated with a hillside 
slope not yet illuminated directly by the early morning 
sun. 


Subsequently the tetroon ascended quite rapidly 
(400 ft. in 2 min. or at a speed of 2 kt.) to a height sufficient 


This sudden ascent took 
place about 2,000 ft. upwind of the ridge itself and as it 


to clear the upcoming ridge. 


could hardly have been thermally induced, must have 
been a purely mechanical consequence of the presence of 
the mountain ridge. 

Figure 7 
traces along flights 14 and 15. 


shows the topography and tetroon height 
The arrows in this figure 
indicate places of maximum (solid arrows) and minimum 
In the case of flight 14, 


which was towed aloft by a radiosonde balloon and flew 


(dashed arrows) tetroon height. 


at an altitude amenable to rather lengthy tracking, the 
tetroon traversed several mountain ridges of varying 
However, the tow release on flight 14 did 
function at the proper altitude with the result that the 


size. not 
tetroon, upon release, slowly sank to its design altitude, 
making the delineation of mountain influences over the 
initial ridge somewhat difficult. On flight 15 there was a 
pronounced tetroon ascent a little upstream from the 
main mountain ridge. On the one hand it would be easy 
On 


the other hand, since the tetroon had previously under- 


to associate this ascent with air flow over the ridge. 


gone a similar rise while over Yucea Flat itself, such an 
association may be erroneous. 

In figure 8 are plotted the height traces (solid lines) 
along flights 14 and 15 as a function of position with 


respect to the nearest mountain ridge. The dotted line 
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shows that in all cases the maximum tetroon height was 
reached upstream from the ridge with the distance up- 
stream being on the order of two times the height of the 
For comparison we have scaled down and 
(dashed lines) 


mountain, 
plotted, in the same figure, 
deduced theoretically by Queney [10] for a mountain of 
The tetroon height traces agree better with 


streamlines 


medium size. 
the low-level Queney streamline than with the high-level 
one. However, it will be apparent to the reader that while 
the tetroons may have a great potential for the study of 
air flow over mountains, these particular flights were not 
Consequently this 


ideally suited for such investigations. 
comparison of observed and theoretical streamlines and 


trajectories will not be carried further. 


CIRCULATIONS IN THE TRANSVERSE 


PLANE 


8. TETROON 


In order to illustrate the nature of tetroon circulations 
in the transverse plane, a straight line was drawn between 
initial and final tetroon position at flight altitude and 
deviations in height and lateral distance from this line 
were computed at minute intervals. Figure 9 shows these 
transverse circulations (looking downstream, or north- 
ward, up the Flat) for the five afternoon flights. Once 
again the tetroon positions are indicated at 5-minute, 
instead of l-minute, intervals for clarity. It is seen that 
basically the transverse tetroon circulations have the same 
vertical and lateral dimensions with an average eddy di- 
ameter in the transverse plane of about 5,000 ft. At Oak 
Ridge, Gifford found eddies of 700-ft. diameter in this 
plane, nearly an order of magnitude smaller than that 
found at Yucea Flat. 

Flights 2, 9, and 10 possess a clockwise circulation in 
the transverse plane whereas flights 1 and 3 


possess a 
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counterclockwise circulation in this plane. Figure 9 shows 
the manner in which the flights were ordered across the 
Flat, from the flight that went farthest to the west (flight 
9) to the flight that went farthest to the east (flight 10). 
It is seen that when the flights are ordered in this way the 
sense of the transverse circulation alternates so that there 
is agreement among the flights as to the position of up- 
ward and downward motion. Thus, although it appears 
unlikely that systematic circulation patterns could be 
maintained for a period of hours in such a region of intense 
circulation, the possibility at least exists that the air 
circulation in Yucea Flat consists of a series of helices. 
perhaps fixed with respect to terrain features, with axes 
running approximately north-south. Figure 10 shows the 
vertical motion pattern which results when an attempt is 
made to draw to the vertical velocity data, obtained from 
the five afternoon tetroon flights, with such helical cireu- 
lation in mind. The axes of the helices, or longitudinal 
rolls, are spaced 15,000-25,000 ft. apart and are positioned 
in such a way that the ascent portion of the helical cireu- 
lation occurs over the western extremity of the Flat but 
apparently not over the eastern extremity. 

It should be emphasized at this point that the vertical 
motion pattern derived from the five afternoon flights, and 
presented in figure 10, is not unique. Along the line of 
the discussion in Section 4, the vertical motion patter 
could also be drawn in the form of transverse rolls oriented 
across Yucca Flat, as in figure 11. The pattern in figur 
11 appears as reasonable as that presented in figure 10 
in fact, most meteorologists, if simply handed the data 
would probably analyze them in the form presented i 
figure 11. The question now arises, which of the circula- 
tion patterns (assuming persistence for either pattern 


presented in figures 10 and 11 is the more basically correct 
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Flight 1 





Flight 3 Flight 10 


Flight 2 
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FIGURE 9. 


positions at 5-min. intervals 


Tetroon circulations in transverse plane (looking downstream or northward up Yucca Fiat) for early-afternoon flights. 
Circulation patterns ordered and spaced across the Flat in agreement with trajectory locatior 


Tetro 
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Ws or can both occur depending upon wind shear and stability There is general agreement in the two traces with indica- 
the conditions? For independent verification we have only tions of height maxima 35,000 and 65,000 ft. north of 
wht fligit 15 which was released during the morning hours Yucca Weather Station and height minima 50,000 and 
0), jusi as the daytime vertical motion systems were becoming 75,000 ft. north thereof. Added support for consideration 
the organized. Figure 12 shows a portion of the height trace of the oscillations as transverse rolls comes from the 
ere of flight 15 (dashed Jine) in comparison with the height approximate agreement in ratio between mixing depth 
up- trace of a hypothetical flight up the middle of the Flat, and (Lagrangian) wavelength obtained from the day- 
“ars as deduced from the vertical motion pattern of figure 11. time tetroon flights (1/3.4) and the ratio of 1/2.7 specified 
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Figure 10.—Isopleths of 5-minute-average vertical motion (knots) Figure 11.—Isopleths of 5-minute-average vertical motion (knots 
Flat derived from the five early-afternoon tetroon in Yueea Flat derived from the five early-afternoon tetroon 
flights assuming transverse-roll circulation patterns 








In Yucca 
flights assuming helical (longitudinal-roll) circulation patterns. 
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Figure 12.—Height trace (solid line) of hypothetical midday flight 
north up Yucea Flat derived from isopleth pattern in figure 11 


2 9, 10) and portion of height trace of 


based on flights 1, 2, 3, 


morning flight 15 (dashed line). 


154] as appropriate to such rolls.* 
pattern of flight 15 


by Scorer [13, p. 
When the 


(fig. 13) is compared with figure 9, however, the agreement 


transverse circulation 
is seen to be very poor because flight 15 has the opposite 
transverse circulation sense to flight 3 and vet the flight 
trajectories almost coincide. This discrepancy is not 
necessarily fatal to the postulation of helical circulation 
patterns because some east-west shift in the axes of the 
With 
the data sample available the problem is moot and 
lies either in the tracking of 
two or more tetroons or in the combination of the Lagran- 


helices might be expected as the sun rises higher. 


solution simultaneous 
gian tetroon data with Eulerian data of the type dis- 
cussed by Richter [11]. 


9. STATISTICS 


While it is not the purpose of this paper to discuss at 
length velocity statistics derived from the Yucea Flat 
tetroon flights, it is useful to consider briefly a few vertical 
indicate their relation to other 
values. Information of this 
It is seen that in the mean for 


velocity statistics and 


velocity component tvpe 
is presented in table 4. 
all the Yueea Flat flights the root mean square vertical 
velocity is 2.4 kt. and the turbulence intensity is 0.25, 
while the ratio of root mean square vertical and longi- 
tudinal, and vertical and transverse velocity components 
is 0.88 and 0.80 respectively. The Wallops Island data, 
involving overwater trajectories, gave values of 0.07, 
Thus, 


desert 


0.80, and 0.46, respectively, for these parameters. 
in comparison with the overwater flights, the 
flights vield «a considerably greater turbulence intensity 
in the vertical and show a lessening of the transverse 
velocity variability with respect to the vertical velocity 
variability, as might be expected. 


*Brought to the attention of the writers by Mr. James E. Caskey, USWB. 
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transverse plane (looking 


for morning flight 15 


Figure 13.—Tetroon circulation § in 


downstream or northward up Yucea Flat 
Tetroon positions at 1-minute intervals. 


There is considerable variation in the above values wit! 
time of day, however. Thus at Yucca Flat the ratio of 
root mean square vertical and longitudinal velocities varies 
from somewhat more than one during midday to about 
one-half at night, whereas the ratio of root mean squar 
vertical and transverse velocities varies considerably less 
This result is in qualitative agreement with the (Eulerian 


findings of Panofsky [9] on a much smaller scale; namely 


that the variability of the transverse velocity is muel 
more dependent upon stability than is that of the longi- 
tudinal velocity so that the ratio of transverse and longi- 
tudinal variability is much larger during the day tha 
at night. 

Correlation coefficients between wind speed and vertica 
motion are also presented in table 4. In the mean fo 
all flights the wind speed is stronger with descending 
motion than with ascending motion so that momentum ls 
being transported downward, as would be anticipated 
On some individual flights, however, momentum is being 
transported upward, particularly during the first day o 
tetroon flights. 

TABLE 4. Vertical velocity statistics derived from Yucca Flat tetroo 


flights. (Height in feet, speed in knots, o is the root mean squar 
R is correlation coefficient.) 


Launched H 
PST 


Flight 


1213 
1346 
1506 
2008 
0635 
1204 
1456 
2036 
O51 


= S~ on ew 


eoonnecsc 


Whom COT Ito 


Average 
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10. CONCLUSION 


While information on the nature and magnitude of 
vertical air motions over Yucca Flat has been obtained 
from these preliminary flights, much still remains to be 
learned about the 3-dimensional air circulation in this 
region. In particular: 

1. It is not apparent whether the vertical circulation 


patterns over Yucca Flat are in the form of transverse or 


longitudinal rolls (helices) or whether both or neither may 
occur in dependence upon wind shear and stability. 

2. The influence of topographic features upon the air 
trajectories has not been clearly delineated. 

3. The transition from upslope to drainage flow, and 
vice versa, has not been successfully explored. 

Future tetroon flights in the Yucca Flat area should 
provide solutions to some of these problems. 
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ABSTRACT 


Further data are presented to indicate the accuracy of the airborne economical radiometer (frequently termed 
radiometersonde when used in a modified radiosonde system) in the measurement of infrared radiation, in view of 


its recent widespread use. Three aspects are discussed. 


economical radiometer with a Suomi ventilated radiometer. 
the net radiation obtained with the economical radiometer in the radiometersonde system. 


The first deals with a nocturnal ground comparison of the 


The second covers an analysis of random errors in 
And, finally, an experi- 


mental in-flight verification of the correctness of the conductivity term in the equations for the economical radiometer 


is discussed. 


1. INSTRUMENT COMPARISON 


A comparison of the Suomi ventilated radiometer [1] 
and the airborne economical radiometer [2, 3] was carried 
out 1.5 meters above homogeneous terrain during July 
1961 at Madison, Wis. 
of both instruments taken at a sample rate of one each 
minute, resulted in 203 readings. The net radiation range 
observed was from 0.0091 to 0.1224 langley per minute. 
Recently Tanner et al. [4] discussed the accuracy of the 


Consecutive nocturnal readings 


economical radiometer in daytime measurements of total 
radiation. 

The statistical results of the frequency distribution of 
the difference in net radiation observed with the airborne 
economical radiometer and the Suomi radiometer are pre- 
sented in figure 1. The class interval for figure 1 is 
.0010 langley per minute. The mean difference between 
net radiation observed with the airborne radiometer and 
the ventilated net radiometer is 0.001059 langley per 
minute with an r.m.s. deviation of 0.001974 langley per 
minute. 

The mean seale factor of the airborne radiometer to the 
ventilated radiometer is 1.01029 with an r.m.s. deviation 
of 0.02017. 
eter readings during nighttime conditions is apparent 
from the results. 


The reproducibility of the airborne radiom- 


2. ERROR ANALYSIS 


In a consideration of the effects of various possible 
random and systematic errors in the measurement of the 
downward- or upward-propagating radiation stream with 
the radiometersonde in balloon ascents, two sources of 
error are suggested by the general equation of Suomi and 
Kuhn {2, 3]. They are random errors in conduction and 
black body radiation resulting from temperature measure- 


lent errors and systematic errors in conduction due to 
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Figure 1.—Frequency distribution of the difference in net radiation 
observed with airborne economical radiometer (radiometersonde) 


and the Suomi ventilated radiometer. 


variations in the vertical dimensions of the instrument. 
We consider, first, errors in the temperature measuring 
system of the top and bottom sensing surfaces of the 
radiometer. 
from equation (1) of Suomi, Staley, and Kuhn [3] may 


The expression for net radiation obtained 


be written: 


Ry ATs—Ts}+2| a {*4 ny} +b] (7 ~7'J 
|e { ett hy a (T.—T}] 
+[e{"3 rn} +4] d\T,—T,) 


T _ TILER 
T.—T,|+Kx di 
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where Py is net radiation in langley per minute; @ is 
0.817107" ly. Ome 
temperature, (°K.); 7, is observed air sensor temperature, 


min. T, is observed top sensor 


(°K.); 7, is observed bottom sensor temperature, (°K.); 
tis time in seconds; A is constant (1.45); and ) is constant 
product of effective specific heat and effective mass/em.? 
of the point-thermistor-Mylar sensor surface experimen- 
tally found from speed of response to be 7 1078 eal./em. 
deg.; a is 0.000005489 langley/minute °C.; 6 is 0.001804 
langley/minute; ¢ is 0.000010979 langlev/minute °C 
d is 0.003608222 langley/minute. 


> and 

The last term on the right side of equation (1) is identi- 

cally zero when the radiometer is tested on the ground. 
Differentiating equation (1) we obtain 

dR +-e)T,+2b 


407’, +- (2a + d\dT, 


40T '—(2a+ce)T,—2b—d\dT,. (2) 


The variance of the net radiation may be expressed as: 


VAR [Ry]J=k? VAR [7,]+h.? VAR [71 


where k, is equal to the first term in brackets on the right 
side of equation (2) and k, is equal to the second bracketed 
T, and 7, are independent of one another, but a 
The 0.2° 


after careful 


term. 
standard error of 0.2° C. is assumed for each. 
C. for the 


checks of the highest accuracy obtainable in reading the 


r.m.s. deviation was obtained 


radiometersonde recorder chart and after extended 
measurements of the normal noise in the recorder system. 
This random error, of course, does not include the fixed 
bias error of the radiosonde receiver-recorder system, but 
the latter effect is of the second order. 

one in which 
70.0° 


the high 


Two cases were solved for equation (3 
the top and bottom radiometer temperatures were 
C. and —40.0° © 
troposphere or low stratosphere, and in the other in which 


respectively, typical of 
these temperatures were 0.0° C. and 15.0° C., typical of 
the lower troposphere. Computation gave a standard 
error of 0.0027 langley per minute in the first case and 
The 
produces an r.m.s. fluctuation of about 1 percent in the 
We 


express the variance as the rate of temperature change in 


0.0041 langley per minute in the second. former 


upper tropospheric net radiation value. can also 
a 50-mb. laver of the atmosphere (see equation (1) In 


Kuhn, Suomi, and Darkow [5]). Since 


OO 1440 g : . 
> | C./day) er ] (4) 


the variance of the left and right sides of this equation 
becomes 
1440 g 


VAR [O07 dt} | ie 


Je VAR (Ry) 


where 07/Odt is the rate of temperature change due to the 
divergence of the net radiation in the vertical, g is the 
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TABLE 1.— Standard error, of net radiation, and of the daily rat of 
radiational temperature change for 50-mb. and 100-mb. atmosph eri¢ 
layers. 


Net radiation oT/or oT /ot 
K-mb. layer 100-mb. lay 


0.0027 ly./min 0.45° C 0.23° C./day 
0.0041 ly./min 0.68° C./day 0.34° C./day 


acceleration of gravity, ¢, is the specific heat of air at 
constant pressure, Ap is the pressure thickness of the 
laver considered, and Ry is the net radiation. 

Solving equation (5) for the two cases cited above 
vields standard errors in the daily rate of radiational 
temperature change for 50-mb. atmospheric layers of 
0.45°C. and 0.68°C., For a 100-mb. layer 
cases I and II give standard errors of 0.23°C. and 0.84°C, 


respectively. 
per day, respectively. These results, of course, involve 
no smoothing, which would tend to reduce the standard 
error. The 50-mb. laver results give a median error of 
0.50°C. for all ascent levels, to be compared with an 
earlier estimate of 0.60°C. [5]. The results described in 
this section are summarized in table 1. 


3. CONDUCTIVITY TERM 


In July 1958, a special balloon ascent was made carrying 
a “dise’”’ type radiometersonde as well as the conventional 
Details of this flight over Madison, 
The importance of the ‘dise” 


radiometersonde. 
Wis. were reported in [6]. 


radiometersonde in this simultaneous ascent becomes 


apparent from figure 2. Since it is exposed to the sam 
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Figure 2.—Cross section of the disc airborne economical radion- 
eter (left) and standard airborne economical radiometer (right 
Cross hatching indicates blackened aluminum sensing surfac 
Ry, is infrared black temperatul 
sensors; subscripted T’s refer to polyethylene, air, top, 4 


radiation current; dots are 
bottom temperatures; upper case P’s to polyethylene ventilati 


shields. 
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E 2.—Comparison of total radiation computed from equation 
for standard radiometersonde and from equation (6) for ‘‘disc’’ 
liometersonde 


ressure (mb.) (stand.) Rot (dise) Riot (disc) — Reos (stand.) 


—0. 035 

0. 030 
—0. 032 
—0. 034 
—0. 035 
—(). 007 
—). 002 
—0. 005 
—0. 001 

0.004 


130 1. 165 
830 0. 800 
797 0. 765 
600 0. 566 
525 0. 490 
440) 0. 433 
340 0. 338 
345 0. 340 
346 0. 345 
347 0. 351 


streams of upward and downward propagating infrared 
radiation, it is clear that in the solution of the radiative 
balance equation for the ‘‘disc’’ radiometersonde there is 
no internal conduction term. Furthermore, in those in- 
stances when the upper and lower polyethylene (air) 
temperatures are equal to the disc temperature there can 
be no conduction between the outer polyethylene venti- 
lation shield and the blackened ‘“disc”’ surface. 

If, then, the total radiation is measured both with the 
“disc” radiometer and the standard radiometersonde, the 
accuracy of the measurement of the internal conduction, 
top polvethylene to top sensing surface and bottom 
polyethylene to bottom sensing surface can be checked. 

The total radiation measured by the ‘‘dise’’ radiometer- 


sonde is given by, 


Rii=20T A —2 [«( mst a )+0 | iT.— T+ oe (6) 


temperature and the con- 
The total radiation 
neasured by the standard radiometersonde is given by, 


Rin=o( T+ 1) —| (A) +4 | tT J 


| e( ae ") +d | [T,—T]+[Knd(T,+-T,)/dt] (7) 


where Ty is the observed ‘‘dise”’ 
stants are defined as in equation (1). 


where all terms are defined as in equation (1). Consider. 
ing only those levels of the ascent where 7,— 74 in equa- 


tion (6) is equal to or less than —1.2°C., we have compu- 
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ted f,., from equations (6) and (7) and compare them in 
table 2. 

The data of the last column of table 2, the deviation of 
the dise from the standard radiometer, display a mean 
value of —0.0117 langley per minute. The r.m.s. devia- 
tion is 0.0206 langley per minute. The mean deviation 
of the total radiation measured by the two instruments 
is approximately 2 percent. Since the same absorptivity 
for the blackened sensing surfaces is used in both instru- 
ments and since the same run polyethylene was used in 
both instruments it is evident that there is no lateral 
conduction loss. For if there were lateral losses the total 
radiation measured by the standard radiometersonde 
would be at best equal to that of the disc for a symmetric 
temperature profile (i.e., the air temperature is the mean 
of the observed top and bottom radiometer temperatures). 
Thus it is believed that the small discrepancy (2 percent) 
in the observed total radiation values springs in part from 
the baseline calibration and in part from the fact that 
T,— Tz in equation (6) was —1.2° C. 

In view of the importance assigned to atmospheric 
measurements of infrared radiation with the radiometer- 
sonde, it is believed that all indications to date suggest 
its continued Studies such as atmospheric flux 
emissivity, tropopause radiation, satellite radiation com- 


use. 


parisons, and others are currently underway. 
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TROPOSPHERIC OZONE: AN AIR POLLUTION PROBLEM 
ARISING IN THE WASHINGTON, D.C. METROPOLITAN AREA 
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D ABSTRACT 

Weather fleck observed on certain varieties of tobacco since 1952 at the U.S. Department of Agriculture Plant 

- Industry Station, Beltsville, Md., may be indicative of an ozone air pollution problem arising in the Washington, D.C 

: area. A study of ozone concentrations and meteorological conditions existing prior to and at the time of the occurrence 
<i of weather fleck indicates that the high ozone concentrations result from the local formation of ozone due to the photo- 
2 chemical reaction of nitrogen dioxide and certain hydrocarbons in the atmosphere of large cities. This reaction is 
J6U) 


possibly enhanced by the prior advection of an enriched-ozone tropospheric air mass 


e it 1. INTRODUCTION 


15 

Ozone has been recognized for some Vears as one of the 
ble more important stratospheric gases. Although its concen- 
a tration is very low compared to concentrations of other 
w atmospheric gases, its physical properties give it an im- 
pris portant role in the heat balance of the stratosphere. 
ork Tropospheric ozone on the other hand was usually con- 
sidered to be unimportant since its concentration (ozone 
. 


density) was only one-fiftieth of that in the stratosphere. 
In terms of concentrations expressed in parts per hundred 
million (pphm) by volume, surface concentrations are 
only 1/300 of those in the stratosphere. This negligible 
surface concentration is due to the continuous destruction 
of ozone near the earth’s surface by oxidation with dust 
and other atmospheric pollutants with which it comes in 
contact 

The source for most of the tropospheric ozone is the 
stratosphere. However in the last decade another source 
of tropospheric ozone has been demonstrated to be of 
importance in atmospheric pollution problems. Labora- 
tory experiments [1] indicate that nitrogen dioxide and 
ivdrocarbons present in certain concentrations when 
radiated by sunlight react to form ozone. The effluent 
‘rom automobile exhausts, industrial plants, oil and gas 
heating, and refuse incinerators injects large amounts of 
hydrocarbons and nitrogen dioxide into the atmosphere. 
It has been found that under the proper meteorological 
and topographical conditions ozone concentrations in 





urban areas may rise to 10 to 20 times the average. The 

Los Angeles Basin “smog” condition is one of the better 

documented examples of the problems which may arise in Ficure | \ typical example of ‘“‘weather fleck”’ injury on the 
ee ee upper surface of a tobacco leaf 
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a large city when it is subjected to an atmosphere contain- 
ing large amounts of ozone and other oxidants. Deteriora- 
tion of materials, damage to crop and other plants, 
reduction of visibility, and deleterious physiological effects 
are greatly increased. 

One of the more sensitive indicators of the onset of the 
ozone-type pollution problem is the appearance of certain 
types of leaf damage to sensitive plants. For example, 
Heggestad and Middleton [2] have shown particular 
varieties of tobacco are subject to a leaf injury known as 
“weather fleck” or “fleck” (fig. 1). The lesions which 
occur on the leaf in the fields were reproduced by ozone 
in chambers with concentrations of about 25 pphm for 
6.5 hours. Ozone at concentrations up to 30 pphm for 
8 hours resulted in fleck and other types of leaf marking 
on 14 of 30 plant species tested [3]. Six of the 14 were 
marked after exposure for 8 hours at less than 10 pphm. 
Since 1952 fleck has been observed on tobacco plants at 
the U.S. Department of Agriculture Plant Industry Sta- 
tion, Beltsville, Md., located about 6 miles northeast of 


Location of stations employed in the ‘rubber-strip’’ network during September and October 1959. 


Washington, D.C. Some fleck may have occurred prio 
to 1952, but there are no records of such occurrence. This 
commencement of fleck damage in this area may indicat 
that the Washington Metropolitan area is being subjecte 
to an increasing number of events of ozone pollutio 
However, because of different topographical and meteor- 
ological conditions it is unlikely that Washington wi 
ever be subjected to smog conditions as severe as thos 
found in Los Angeles. 

The examination of the meteorological conditions exist 
ing at the time of high ozone concentration was begun! 
1958 as a joint project of the U.S. Weather Bureau and th 
Department of Agriculture Plant Industry Station. Thi 
results of that project have been reported [4]. The proje 
was continued in 1959 and a small network of ozone-mea* 
uring stations surrounding Washington was added. T! 
results for the 1959 period are the subject of this paper. 


2. RUBBER-STRIP NETWORK 


In the 1959 tobacco growing season a small network 
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statio s was established at which rubber strips were ex- 
posed daily. 
mont!.s of September and October. 


These stations were operated during the 
The locations of the 
statios with respect to the district of Columbia are shown 
in figure 2. 

The use of rubber strips as instruments for the detection 
of ozone, described previously [5, 6], will be discussed here 
only briefly. When rubber is under tension, ozone causes 
a particular type of cracking (fig. 3). Laboratory investi- 
gations have shown that this cracking is due only to ozone. 
However, because of the number of parameters which in- 
fluence the action of ozone on rubber, the determination 
of atmospheric ozone is more qualitative than quantitative. 
Parameters which are important are: (1) Nature of formu- 
lation and degree of stress of the rubber compound. (2) 
Concentration of ozone. (3) Time and method of expo- 


sure, including shading. (4) Amount of ventilation re- 


ceived by rubber. (5) Temperature and humidity. 

Some subjectivity is encountered in determining the de- 
eree of cracking as the strips are examined under a low- 
power microscope. The results depend somewhat on the 
visual acuity of the observer. 

Rubber strips, 40 mm. S min. 1 mm., were bent and 
placed under stress by mounting them in spring clips which 
at most stations were suspended within inverted plastic- 
coated paper drinking cups. The bottoms of the cups were 
cut half open and interiorly baffled to provide both shade 
and natural ventilation. Exposures at Centreville from 
October 15 were made inside a standard weather instru- 
ment shelter. Two rubber strips were exposed side by side 
nthe same clip at Brandywine and at the Weather Bureau 
Central Office ; only one strip was exposed at other stations. 
At Beltsville the spring clips were suspended below a hori- 
zontal plank and at Upper Marlboro two strips, one below 
a standard shelter, the other inside the shelter, were ex- 
posed. The depth of individual cracks was measured 
under low-power magnification and totaled for two freshly 
cut edges of each strip by the procedure emploved by 
Haagen-Smit and Bradley [5]. Where two strips were ex- 
posed, the average total was obtained. 

The rubber-strip network operated from September 10 
through October 25, 1959. All strips from Beltsville, Cen- 
ireville, and Upper Marlboro were evaluated for total 
crack depth. 
the entire network were evaluated for only September 22 
to October 7, 1959. 


Because of cost considerations, strips from 


In addition to the rubber strips, the U.S. Department 
of Agriculture Plant Industry Station operated an auto- 
matic chemical ozone meter similar in principle to the one 
ised in 1958 [2, 4]. Because of uncertainties regarding the 
absolute values of ozone concentration indicated by the 
recorder, a relative scale of 0 to 100 is employed in this 


paper to express the daily peak concentration. During 


the 1959 season, the ozone meter operated satisfactorily 
irom September 21 to October 7, 
and 6. 


except on October 3 
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Figure 3.—Typical cracking of a rubber strip under tension when 


t xposed to atmospheric ozone, 


3. HIGH-OZONE DAYS 


In the paper already cited [4] a high-ozone day is defined 
as one on which an 8-hour average of ozone concentration 
was 31 or greater on our relative scale of 0 to 100. On this 
scale the average daily peak value was 29 and the median 
was 22, whereas on high-ozone days the peaks ranged from 
62 to 100. 
24, with relative peak concentrations of 62 and 68 at 1110 


Two days in the 1959 series, September 23 and 


and 1215 Est, respectively, qualified as high-ozone days 
by our definition. On the remaining 15 days, the daily 


peak values ranged from 14 to 46. 


4. CORRELATION OF RUBBER CRACKING WITH 
AVERAGE OZONE CONCENTRATION 


On high-ozone days in 1958, the times of peak values 
ranged from 0955 to 1505 est. When the peak occurred 
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FIGURE 4.—Scattergram of daily peak and average ozone values 
0900-1700 EsT) on a relative scale at Beltsville, Md. 
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early, the rise from near-zero values at night was abrupt, 
compared with the rate of falling off in the late afternoon. 
It is known that the cracking of rubber depends upon 
both the concentration of ozone and the rate of air move- 





ey 
°o 
r 


ment past the rubber strip. Since nighttime values of 


ozone at Beltsville are nearly zero except on windy nights 
°9-23 


a ee 


and wind speeds tend to have a maximum soon after 


ow 
°o 


midday, it was supposed that the amount of rubber 

— . . . . . ' 
cracking in a strip exposed for 24 hours was a function °9-25 
of its exposure to ozone and wind speed during the 8-hour 


period between 0900 and 1700. The correlation between 


nN 
° 





daily peak values and average values for this time span 


is fairly good (fig. 4). On October 4 and 5 the peak 


values occurred as relatively brief spikes, consistent with 


ray 


the high ratio of peak to average values. The ratios 
of peak to average 11-hour concentration (0900 to 2000) 
and to average 24-hour concentration (midnight to mid- 


night) were found in approximate ratios 2:1 and 7:2, 


DAILY AVERAGE OZONE LEVEL, 0900-1700 E.S T. (RELATIVE SCALE) 
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tion at Beltsville are compared directly, little useful 


respectively, with somewhat more scatter. 


correlation is observed (fig. 5). The cluster of points at : ee wer ae ’ 
‘ = heray é I FIGURE 6.—Seattergram of ratio of total crack depth to averae 
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wind s eed (0900-1700) are compared with average ozone 


levels, « reasonable correlation is found (fig. 6). 


5. WIND SPEED CORRECTION 


The wind speed used to correct for natural ventilation 
of the exposed rubber strips is the average of the speeds 
measured at Chantilly, Va. (site of Dulles International 
Airport) and Upper Marlboro, Md. (fig. 2). The anemom- 
eter heights are 20 and 30 feet, respectively. 
measured at the Washington National Airport were not 
included because of a change in anemometer height from 
The allowance for ven- 


The speeds 


” to 115 feet on September 30. 
tilation made in this paper is of course an approximation. 
Hence significance must not be given to small differences 
We 


shall consider only the highest ratios and the larger differ- 


n the daily ratios of crack depth to wind speed. 


ences in the subsequent analysis. 
6. COMPARISON OF 1958 AND 1959 RESULTS 


We consider first the two high-ozone days September 
23 and 24 and compare them with the 5 high-ozone days 
11958. The correspondence of the ensemble of weather 
parameters reported for the 1958 cases with the values of 
the weather parameters on these two daies is generally 
worable. (See table 1). Light winds from the direc- 
m of the metropolis together with high solar irradiation 

‘similar to 1958. Lower humidities in the tobacco 
eld and higher temperatures at both tobacco field and 
irport are observed. 

Severe weather fleck was observed at Beltsville on two 
Physiologica! 


lays in 1959, July 27 and September 25. 


Some meteorological parameters on high-ozone days, 1958 


and 1959, Reltsville, Md. 
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{at Washington National Airport, 14 mi. south-southwest of Beltsville. 
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Figure 7.— Areal 
depth for the metropolitan and suburban areas of Washington, 
D.C. on September 24, 1960. 


distribution of wind direction and total crack 


insult by ozone is presumed to have occurred the day 


before. The later date corresponds to the date of the 
highest peak ozone concentration in the period of recorder 
operations, September 24. The time of the peak value 
on July 26 is unknown, but if we take 1530, the corres- 


The 
differences from the 1958 experience on July 26 include 


pondence of the weather ensemble is good (table 1) 


higher air temperature, smaller temperature range, and 
the occurence of slight precipitation on the preceding 
day, none of which represents basic changes. The com- 
bined 1958 and 1959 results presented in table 1 should 
therefore form a basis for further analyses of high-ozone 
days near a metropolis and assist in the more difficult 
analyses of other days with less ozone but still higher ozone 
content than normal. 


7. SOURCE OF OZONE 


An examination of available ozone data and meteoro- 
logical charts was made to ascertain the source of the high 
concentration of ozone observed on September 23 and 24 
at Beltsville, Md. 


city air with the proper concentration of nitrogen dioxide 


If ozone is produced by irradiation of 


and hydrocarbons, then on days when light winds cross 
the city steadily and solar irradiation is high, we would 
expect an increase in rubber cracking on the downwind 
side of the city. 

On September 24, a day of appropriate wind and irradi- 
ation, the results indicate that this was indeed the case 


(fig. 7). Most frequently reported surface wind direction 
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TABLE 2 
Date 


Central Beltsville Centreville 


September 1959 
office 


and extremes during the period 0900 to 1700 Est, Sep- 
tember 24, together with wind directions reported for 300 
meters above sea level at Silver Hill, show a consistent 
flow from south or southwest across the metropolitan 
Flow from this octant will usually exclude Balti- 
more as a source region. Rubber-cracking ratio at Belts- 
ville, downwind of the city, was 2.29, compared with 1.04, 
1.01, and 0.98 upwind or crosswind at Herndon, Centre- 
ville, and Brandywine, respectively (see table 2). The 
Central Office and Baltimore strips were exposed 30 inches 
above roof tops, where wind speeds are known to be 
higher than at 30 inches above the ground. Hence the 
values 1.24 and 1.87 are probably both high. The rela- 
tively low value at the Central Office within the city might 
be ascribed to insufficient irradiation of the ozone pre- 
cursor material in its travel from the outskirts of the 


urea. 


metropolitan area to the south toward the site or to the 
presence of ozone sinks inside the city. A similar but 
less-pronounced gradient of rubber-cracking ratios was 
observed on September 23, the other high-ozone day. 

An examination of the areal average of the rubber- 
cracking ratios (table 2) indicates that the background 
level of ozone began to rise significantly on September 
20. A peak value of the ratio was reached on September 
25, and a steady decrease until the end of the month 
followed. The steady rise in the ozone level may have 
been due to the advection of ozone-enriched tropospheric 
air over the area. 

The source of the tropospheric air arriving over Wash- 
ington during September 20 to 24 was determined from 
air trajectories on the 850-mb. surface. From the 20th 
to the 22d the trajectories indicated a polar source region 
in the vicinity of Hudson Bay, while the period from the 
23d to the 24th indicated a source region over the eastern 
portion of the Gulf of Mexico. An inspection of the 
maximum wind charts indicated that while each of the 
air Masses was in these source regions a jet stream axis 
was also present over the region. Unfortunately, meas- 
urements of the vertical distribution of ozone for these 
areas are not available for one to determine whether a 
downward transport of ozone from the stratosphere to 


Herndon 


Individual station and areal ave rage values of the ratio of total crack de pth to wind speed Jor Se ple mber 15-30, 19 


Ratio of total crack depth to average wind speed 


Davidson Brar‘ywine | Silver Upper 


Baltimore 
ville Marlboro 


the troposphere occurred during these 


previous experience {7] has indicated that such processes 


are possible near the jet stream axis. Investigations 0 
the vertical velocities in the vicinity of jet streams [84 
have shown that they may serve as a mechanism fo 
carrying out vertical exchange processes in this regio: 
Endlich [8] found that the strongest vertical speeds wer 
located near the jet axis. Vertical speeds of 5 em. see 
were seldom found farther than about 300 miles from th 
axis. Endlich and McLean [9], in a detailed analysis, 
the jet stream core based upon aircraft flights into tl 
jet stream, reported that downward vertical speeds 


‘ times encountered in t! 


excess of 1 m. sec.~' were at 
isotach centers of the jet. 


Measurements of the vertical distribution of ozo: 


in the vicinity of the jet axis were reported by Morela 


[7]. The measurements were obtained during a serie: 


Denve 
An example of an ascension ma 


of intercomparison ozonesonde flights made at 
Colo., in April 1959. 
in the vicinity of the jet stream is shown in figure » 


The axis of the jet, 4 hours before release time, wa 


located about 100 miles northwest of Denver (fig. 9 
The aerological data (fig. 8) taken simultaneously wit 
the ozonesonde ascent denoted a well-defined tropopaus 
and a maximum wind speed of 40 m. sec.~! (78 kt. 

11 km. The ozone concentration (ozone density) bega 
« marked increase at about 9 km. indicating that t! 
ozone had been injected into the troposphere from th 
stratosphere. In contrast, the vertical distribution ¢ 
ozone shown in figure 10 was obtained when the je 
stream was about 400 miles east of Denver (fig. 1! 
Under these conditions the ozone concentration did no 
begin its increase until 4 km. above the tropopaus 
Therefore preliminary evidence suggests that the I 


jection of stratospheric ozone may occur in the vicinit! 


of the jet stream axis. 
8. CONCLUSIONS 
Thus from all available evidence, which admitted! 


is Meager, it appears that the occurrence of high 020! 
concentration in the Washington area may be the rest 
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(1) the increase of the general back- 


of two 
ground level of ozone concentration by the advection 


processes - 


of an enriched-ozone tropospheric air mass; (2) the local 
formation of ozone under stable atmospheric conditions 
at low levels due to the irradiation of precursor material 
as described previously. Evidence of an ozone gradient 
across the city with high concentration of surface ozone 
enhanced downwind supports the hypothesis that the 
metropolis is the dominant local source of ozone formation. 
The may 
sufficiently 
plant dumage observed at Beltsville [2]. 


combination of these two processes cause 
high ozone concentration to produce the 
Whether one 
of these processes acting alone might be sufficient to 
cause the fleck damzege cannot be determined with the 
umount of date now available. However, if more ozone 
observations become available, a more detailed analysis 
of the 


plant damage and air-pollution potential can be under- 


ozone-meteorological conditions conducive to 


tuken. 
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Over a half century after the discovery of its existence, 
the tropopause is still one of the challenging mysteries to 
indicate that 
the vertical temperature gradient of the air often changes 


meteorologists. Radiosonde observations 
rather abruptly at the tropopause, although the apparent 
sharp break of the temperature gradient may be artificial 
because of the fact that straight line segments between 
reported temperature observations are employed to repre- 
sent the true, and perhaps smoother, temperature dis- 
Even if the kink of the temperature curve at 
ihe tropopause is removed, the fact remains that the 


tribution. 


temperature gradient changes rapidly within a rather 
narrow zone immediately above and below the point desig- 
nated as the tropopause. The mystery about the tropo- 
pause is, therefore, not primarily its location, but rather 
the sharpness of the change of the temperature gradient 
that is associated with it. The radiative processes and 
the other large-scale heat transport processes that deter- 
mine the mean atmospheric temperature distribution 
probably also determine the large-scale features of the 
irepopause, such as its general location and its seasonal 
id large-scale spatial variations. But what is the reason 
for the often observed sharp change of the temperature 
gradient at What are the physical 
processes that are responsible for its creation, maintenance, 


the tropopause? 


transformation, and destruction? 
Few 


suggestions concerning the creation of a sharp tropopause 


There are as vet no answers to these questions. 


r its maintenance have been put forward (i.e. Moller 
j], Staley [8], Kaplan [3]) and they are hypothetical in 
that they are based on other assumptions. 

Still more basic and serious than this, it seems to this 
writer, is the fact that we do not even have an adequate 
bservational study of the tropopause. We have only a 
gross and simplified picture of the tropopause and we do 
iot know how a particular tropopause evolves with time. 

There are several reasons why ordinary radiosonde 
observations fail to provide a detailed picture of the 
‘ropopause. First of all, the ordinary practice of report- 
ug alternatively the humidity and temperature prohibits 
‘continuous description of temperature. 
result of the combined effect of the large ascensional rate 


Second, as a 


ithe balloon and the time lag of the sensing thermistor, 
‘he reported temperature actually represents an average 
‘mperature over a rather large distance in the vertical 
‘low the location of the reported temperature.? In 
dition, the representation of the temperature trans- 
‘ssion in each baroswitch contact by a short trace on the 
‘cording paper often makes it undesirable to obtain more 


than one or two temperature readings in each baroswitch 
At the 200-mb. level, which is about the level 
of the tropopause at middle latitudes, each contact repre- 
sents a vertical distance of about 450 feet. Therefore, 
the maximum meaningful information that one may derive 
from the ordinary radiosonde is about one or two tempera- 


contact. 


ture readings for every 450 feet at the tropopause laver. 
This is not enough to delineate the detailed structure of 
the tropopause. 

Equally as important as information about the detailed 
structure of the tropopause at one particular moment, 
is knowledge of the evolution of the same tropopause with 
time. Efforts to follow the life history of a tropopause 
have met with little success due to the fact that the ob- 
servations now available at fixed locations are not meant 
and are not suitable for this kind of study. Only by fol- 
lowing and making observations of one and the same tropo- 
pause, can the life history of the tropopause be ascertained. 

In view of the above inadequacies of the present radio- 
sonde system, it is hereby suggested that instrumented 
constant level balloons be employed to follow and to probe 
the tropopause. This suggestion was prompted by the 
following two observations: 

1. At both low and high latitudes the tropopause is 
often well defined and coincides with the minimum tem- 
perature. It seems, therefore, that the tropopause is a 
continuous surface of minimum temperature with con- 
siderable horizontal extent at those latitudes. 

2. In the last decade or so, constant level balloons have 
been successfully used to probe the air flow of upper 
troposphere or lower stratosphere (Moore, Smith, and 
Gaalswvk [6], Neiburger and Angell {7)}). 
were set to float along a certain isobaric surface. 


The balloons 


If the tropopause is a surface of minimum temperature, 
then in principle it should be possible to set the balloons 
to stay and to float at this surface of minimum tem- 
perature. This can be done, for example, by having a 
series of simultaneous temperature observations (3 or 
more) from thermistors attached an equal distance apart 
on a line dangling from the balloon, and by having the 
floating level of the balloon regulated by the requirement 
that the middle thermistor should register the minimum 
If due to the ascending or descending of 


temperature.* 


This correspondence is an abridged and revised version of a scientific report [1] of the 
same title written under Contract A F19(604)-6146 sponsored by Geophysics Research 
Directorate, Air Research and Development Command 

? For one type of U.S. radiosonde, the vertical distance of averaging is estimated to be 
no less than 350 feet (Chiu [1] 

‘ Since the balloon package is floating with the air, an aspirator would be necessary to 
keep each of the temperature sensing elements ventilated in order to prevent radiation 
errors 
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the balloons this requirement is not met, then the com- 
parison between the temperature registered by the middle 
thermistor and its immediate neighboring thermistors 
will set a mechanism into operation (e.g., the dropping 
of the ballast or the releasing of the gas from the balloon) to 
bring the balloon up or down to the desired level. In 
this manner it should be possible to control the balloon, 
causing it to float and to oscillate about the tropopause 
with a certain preset limit (this limit is determined by the 
spacing between the thermistors and by the sensitivity 
of the height control mechanism). 

The above described observations, if successful, could 
provide us with a detailed temperature structure of the 
tropopause laver (a tropopause layer is here interpreted 
us meaning a layer of certain thickness in which the tropo- 
There are many other useful pieces of 
For example, 


pause is found).* 
information which we should like to have. 
if the humidity can be measured with certain desired ac- 
curacy then we may replace all the thermistors in the above 
experiment by humidity elements and get a detailed struc- 
ture of the vertical humidity distribution. The im- 
portance of this kind of information is obvious. The 
difficulty is, of course, whether equipment can be found 
which will measure the humidity with adequate accuracy 
at the tropopause level. 

Another very useful piece of information that might be 
obtained from the floating balloons is the heating or cooling 
of air in the tropopause laver by radiative processes. 
This may be done, for example, by installing three of 
Suomi’s [9] net radiometers at, above, and below the tropo- 
pause. From these measurements the divergence of radi- 
ative flux, and therefore the temperature change, above 
and below the tropopause could be easily calculated. In 
this way it might be possible to see just how the radiative 
processes are affecting the detailed structure of the tropo- 
pause, providing that the balloons are traveling with the 
same alr parcel, 

It is also possible that the radiation measurements at 
the tropopause may be used to settle one of the very 


important questions about the tropopause; i.e., is the 


tropopause a material surface? ® By material surface is 
meant one consisting of the same air particles as time 
progresses. If the tropopause is a material surface and 
the ballon follows the horizontal motion of the air, then 
the balloon will stay in the same parcel of air as it travels. 
In that case there should be no advection of heat into the 
The change of the poten- 

the 
then be mainly due to radiative processes. 


temperature sensing elements. 
thermistors would 


By comparing 


tial temperature registered by 


the mezsured potential temperature changes from the 
radiometers, on the one hand, and those from the thermis- 
tors, on the other, it should be possible then to ascertain 
whether or not the tropopause is a material surface. 

* According to the people 
Bedford, Mass., 
permits up to 25 temperature readings on a 2,000-foot line 

In the 


been expressed 


it Instrument Research Laboratory, Geophysics Research 


Directorate an encoder-thermistor arrangement is now available that 
iterature, opinions both for and against the idea of a material tropopause have 


Lowell [4], Danielsen [2 
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An alternate scheme of carrying out the observation js 
to set the balloon to float freely near the tropopause. — This 
scheme was suggested to the writer by the Air Foree 
balloon design and operation groups at Geophysics Re- 
search Directorate, Bedford, Mass., when they were 
consulted on the balloon flight problem. 
three people [10], ‘the balloon could support a line about 
2,000 feet long beneath it, 
appropriate intervals along the line. 
be set to float at or just below the expected height of the 
tropopause, and its altitude could be adjusted upward in 
small increments by the dropping of ballast.’’ In this 
way, there is a fair chance that the line dangling from the 


According to 


with thermistors spaced at 
The balloon could 


balloon will traverse the tropopause. 

At this stage, the suggested new technique of observing 
the tropopause is only an idea. Funds for its experimen- 
tation must be acquired. Much effort then 
made to design and develop the balloon and its instru- 
Only after many 
But in 


view of our meager knowledge about the tropopause, if 


must be 


mentation and to test them in flights 
flight tests can one begin to assess its usefulness. 


this correspondence can serve in any way to stimulate 
interest and effort into the research of the tropopause 
then it already has served a very useful purpose. 
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THE WEATHER AND CIRCULATION OF MAY 1961 
Persistent Cool Weather in the United States 


ROBERT H. GELHARD 


Extended Fore 


1. HIGHLIGHTS 


The outstanding feature of the weather of May 1961 
was the persistence of the cold regime of the preceding 
month [1] throughout most of the United States east of 
the Rocky Mountains and also in the Pacific Northwest. 
In the East stations from St. Mo., 
eastward through Ohio Valley to Youngstown, 
Ohio reported the coldest average temperatures on record 
for May. 


tures retarded the maturity of crops through most of 


numerous 
the 


Louis, 


In the Far West unseasonably cold tempera- 


the Great Valley of California and below normal tem- 
peratures extended into extreme southern California. 
A major contribution to the low mean temperatures 
for May came from an extremely cold polar air mass 
that blanketed the eastern half of the country during 
the 
record low temperatures were recorded throughout the 
Middle West 


and 


the last week of month. Beginning on the 26th, 
and eastward, with frost as far south as 
Ohio Michigan. Daily 
temperature anomalies ranged from more than 20° F. 
below normal in New England and the Mid-West to 
10° F 


period. 


Kentucky snow in and 


below normal through the Gulf States during this 


Extremely heavy precipitation from the Central Plains 
eastward through the Ohio Valley resulted in major floods 
throughout the month in most of the Arkansas, Missouri, 
Ohio, and Middle Mississippi River watersheds where 
amounts ranged from more than 8 inches at Concordia, 
On 
the other hand, drought conditions that prevailed in 


Kans. to a record 12.22 inches at Evansville, Ind. 


April 1961 through the Southwest from Texas to southern 


California continued during May except in extreme 


western Texas. There was no measurable precipitation 
in the extreme Southwest or at Corpus Christi, Tex., 


inew record for May at that station. 


2. THE GENERAL CIRCULATION 


The general circulation for May 1961 at 700 mb. 


ig. 1) does not show many unusual features in spite of 


the extremely anomalous weather observed over most 
of the United States. 
easter 


\tlant 


Blocking, which had dominated 
Canada during April [1] and the northeastern 
during March [2], retrograded rapidly to the 


Alaskan area where mean height anomalies increased 
210 ft. from April (fig. 2). 
anomalies decreased by 570 ft. near Baffin Island, repre- 


At the same time height 


senting a rather abrupt return to normal in that area. 
In the Southwest Pacific filling of the deep Asiatic 
coastal trough observed in April and intensification of 
the subtropical High downstream were predominantly 
climatological as shown by the small anomalous changes 
(fig. 2). However, this normal change, concurrent with 
deepening east of Kamchatka and blocking to the north, 
produced a zone of confluence south of Kamchatka which 
was the dominant feature of the circulation. Figure 3 
that the 700-mb. 


this region were from 4° to 5° of latitude north of their 


shows maximum mean westerlies in 
normal position (A) and 9 meters per second faster than 
normal (B). 


this high-index Pacific circulation are shown by a diminu- 


Downstream, the effects of blocking on 


tion of wind speeds in the broad difluent area of western 
North America and by the southward displacement (from 
normal) of the 700-mb. jet stream through the extreme 
southwestern United States. 

Elsewhere the anomalies of the general circulation were 
160° ft. 
month-to-month changes (fig. 2) being indicative of a 
the 


within of the long-period normal, the largest 


return to normal from more anomalous state of 


April [1]. 
3. TEMPERATURE 


The cool temperature regime of April continued through 
May with nearly all of the contiguous United States aver- 
aging below normal except Texas and the Rocky Moun- 
tain States (fig. 4). In much of the Mid-West and Ohio 
Valley the persistent low temperatures of April and May 
1961 combined to make this period the coldest since 1907. 
The cool regime in the Pacific States and the Great Basin, 
while not of record intensity, was extremely persistent 
particularly in the Great Valley of southern California. 
The number of days in the month during which the aver- 
age daily temperature exceeded the normal ranged from 
only six at Sacramento to four at San Diego. Santa 
Maria recorded only one day above normal and that by 
only 1° F. 

The unusual month-to-month persistence of temperature 
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Ficure |. 


Mean 700-imnb. contours (solid) and height departures from normal (dotted) 


both in tens of feet) for May 1961 Above norma 


heights in Canada and Alaska were associated with frequent anticyclogenesis and cool weather east of the Rocky Mountains 


is well illustrated by figure 5, the change in temperature 
anomaly class (out of five classes) from April to May. 
In a distribution of 100 cities, 37 remained in exactly the 
same class as the previous month, covering the area shown 
No city changed by as many as 3 or 
The 
latter figure compares with 59 expected by chance, 32 
observed from March to April 1961 [1], and 62 for the 
average April-May change from 1942 to 1957 [4]. 
Although the pattern of temperature anomaly follows 


by crosshatching. 
4 classes, and 79 changed by no more than 1 class. 


closely the negative height anomaly at 700 mb. the rela- 


tive magnitude of the temperature anomaly was mue! 
greater than that of the height anomaly. Examination 0! 
the 15-day components of the monthly 700-mb. height 
mean, figure 6, shows a two-fold effect operating to mail 
During the first 


tain cool temperatures in the Far West. 
wit! 


half of the month low temperatures were associated 
anomalous northerly flow from a blocking High centered 
over Alaska (fig. 6A), and with the main westerlics dis 
placed well south of 40° N. During the latter half of th 
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Figurk 2.— Change in 700-mb. height anomalies (tens of feet) from 


April 1961 to May 1961. 
ada is represented by the largest change area. 


Relaxation of blocking in eastern Can- 


month temperatures modified considerably in the North- 
west (fig. 8 C, E) as the blocking over Alaska relaxed (fig. 
7), and the circulation returned to the dominant spring 
pattern of troughs near the east and west coasts and a 
strong full-latitude ridge through mid-continent (fig. 6B). 
The extreme negative tilt (northwest to southeast) of the 
west coast trough is commonly associated with the south- 
ward displacement of mid-latitude westerlies and helped 
maintain cool temperatures in southern California. 

In contrast to the persistent temperature regime of the 
Far West, temperatures east of the Continental Divide, 
showed 


normal, more 


This variability can be related to 


although predominantly below 
variability (fig. 8). 
the differences in the 700-mb. bimonthly circulations 
hig. 7 

During the first half of the month most of the United 
States was dominated by a broad cyclonic circulation aug- 
nented by strong northerly flow in western Canada asso- 
cated with a blocking High over Alaska. Following an 
nitial outbreak of cold air which affected most of the 
ountry except the Gulf States, succeeding polar air 
lasses followed a more northerly track with a warming 
irend east of the trough in central United States. The 
track of the “glancing Highs” associated with these shallow 
polar outbreaks is shown in figure 9A by the arrows 
through the Great Lakes and eastern Canada. Examina- 
tion of the 5-day mean charts shows a transition period 
tear mid-month during which the maximum-curvature 
‘rough in central Canada joined with the shallow trough 
i the Northern Plains (fig. 6A) to form a full-latitude 
‘rough near the east coast, and the blocking High over 
Alaska moved east to form a full-latitude ridge just east 
of the Continental Divide (fig. 6B). Associated with the 
evelopment of this large-amplitude trough and ridge 


‘stem was rapid warming throughout the mountain 
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Figure 3.—(A) Mean 700-mb. isotachs and (B 
normal wind speed (both in meters per second 


departure 
for May 
Solid arrows in (A) indicate primary axes of the westerly jet 
stream, and dashed arrows indicate their normal May positions 
The westerlies were displaced well south of the normal position 
over western United States, and an important secondary maxi- 


mum was located over northeastern Canada 


area of western North America (fig. 8G). A succession of 
cold polar anticyclones moved into the eastern half of the 
country, as shown by the southernmost track in figure 9A, 
resulting in frost and light snow in the Mid-West and 


record low minimum temperatures throughout the East 


4. PRECIPITATION 


Of particular interest is the unusually heavy precipita- 
the Central Plains and Middle 
Mississippi watershed and resulted in floods throughout 


tion which occurred in 


the area most of the month (fig. 10). Heaviest rains fell 


during the first half of the month (fig. 8 B, D) associated 
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Figure 4 Departure of average surface temperature ( 
normal for May 1961. (From [3} 
normal (shaded) in the 


in the remainder of the 
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Temperature averaged above 
Mountain area but below 





tock | Figure 5.— The number of classes (out of five) the surface tempera- 
‘ky norm: : 

eae citi ture anomaly changed from April to May 1961 

nation. 


The hatched area 
represents the area of no change 























MAY 16-30 minus 1-15. 1961 








FIGURE 7. 


Change in 15-day mean 700-mb. height anomaly fro 


May 1-15 to May 16 


30, 1961, in tens of feet. Large intra- 
monthly changes occurred over most of North America except t! 
west coast of United States 


FIGURE 6. 
MAY 16-30, 196! 





Fifteen-day mean 700-mb. height contours an 
tures from normal (both in tens of feet 
and (B) May 16-30, 1961. 
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depar- 
for (A) May 1-1 196! 
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Deporture of Average Temperature trom Normal (°F) tor the Week Ending 


Midmght ist May 14, 196 
































Total Precipitation, inches, for the Week Ending 


Midmght Ist. May 14, 196! 












































Figure 8 Departure of average surface temperaure from normal (° F.) and total precipitation (inches) for the weekly periods ending 
(A,B) May 7, (C,D) May 14, (E,F) May 21, and (G,H) May 28, 1961. (From [6] 
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B. CYCLONES MAY 196] 











Fiat RE 4 Frequeney of (A 


squares at 45° N.) 


anticyelone and (B) cyclone passages 
during May 1961. 
are indicated by solid arrows, secondary tracks by dashed arrows 


within 5 Primary tracks 


Frequent weak Lows in the Northern Plains were of secondary 
importance compared to major storm development in Central 
Plains 


with the shallow trough through the Plains States and the 
The effect of the 
depressed westerlies in restricting the northward intrusion 
of Gulf moisture is shown indirectly in the frequency and 
prevailing tracks of daily cyclones for the month (fig. 9B). 


depressed axis of prevailing westerlies. 


Although a maximum number of closed Lows was observed 
in the northwestern Plains, few developed into well 
organized cyclonic systems, while the major developments 
occurred just north of the prevailing westerlies in the 
Central Plains. 
limit of Gulf moisture was the range of precipitation 


An interesting effect of the northward 


through Indiana from a record 12.22 inches at Evansville 
in the south to only 1.19 inches at South Bend in the north, 


the dryest there since 1939. Such sharp gradients are 





Percentage of Normal Precipstohon 
May, 196! 


























Figure 10.—Percentage of normal percipitation for May 196! 


(From [3]. No measurable precipitation in parts of the South- 


west contrasted with above normal amounts in many other areas 


frequently associated with marked confluence of the mean 
700-mb. contours, such as prevailed over Indiana this 
month (fig. 1). Precipitation was normal to above normal 
in the Pacific Northwest under the influence of the weak 
coastal trough despite the weakness (relative to normal 
of the westerlies at mid-latitudes. 

During the week ending May 15, a major storm moved 
out of the lower Mississippi Valley and produced moderate 
to heavy precipitation over most of the eastern half of the 
country. Intense deepening of the parent storm as it 
moved into eastern Canada resulted in a fresh polar out- 
break of cold air reaching to the Gulf States, and a rather 
intense wave developed in response to the sharp tempera- 
ture gradient near northern Florida. However, a sudden 
northward shift of the jet stream in response to the deep 
cvelonic vortex near Hudson Bay and anticvelogenesis to 
the south over the Great Lakes resulted in cutting off the 
wave near Florida from the main westerlies. The major 
portion of the month’s precipitation in the coastal area of 
South Carolina to Virginia was produced by the slow 
northward migration of this cut-off coastal Low from 
May 10 to May 13. 

Heaviest precipitation in the Northeast occurred late 
in the month in conjunction with wave development as a 
cold front moved off the Middle Atlantic coast in advane 
of the extremely cold polar air mass mentioned earlier 
Subsequent baroclinic deepening in the mean east coast 
trough resulted in precipitation amounts in excess of 4 
inches in parts of New England (fig. 8H). 

As is often the case, precipitation, or the lack of it, 's 
difficult to specify. In the Southwest it seems likely that 
the strong low-latitude westerlies operated to produce & 
rain shadow effect and to prevent moisture from being 
introduced by large-scale southeasterly components of the 
More difficult to explain is the con- 
tinued drought in eastern Texas. Here lack of moisture 
at high levels and anticyclonic vorticity just south of the 


general circulation. 


jet stream limited precipitation to scattered showers 








Aves 1961 


During the latter half of the 
a strong easterly flow at low levels associated with 


throuch most of the period. 
mont 
, series of large polar anticyclones moving along the 
southern track of figure 9A resulted in a considerable 
mount of up-slope precipitation in western Texas. An 
nteresting note on the discontinuous nature of precipita- 
tion is the difference between the monthly totals at Hous- 
ion, Tex. of 0.63 inches at the city office and 3.59 inches 
it the Airport Station. 

In an effort to specify more clearly the month’s precipi- 
tation, a study was made of the 600-mb. vertical motion 
harts as computed twice daily from the baroclinic model 
currently used by the Numerical Weather Prediction 
Branch (NWP). Values were read at 5° intervals of 
latitude and longitude, and charts were constructed of 
veekly mean values and mean positive values only for 
periods corresponding to those in figure 8, and mean 
values and mean positive values only for the period 
orresponding to figure 10. The results show little corre- 
spondence between the observed precipitation patterns 
and the patterns of mean values of vertical motion on 
ither a weekly or a monthly basis. On the other hand, 


the patterns of positive values only show a fair corre- 
spondence with the observed precipitation in some areas. 
the dry 


reas of Texas and the Great Lakes show definite maxima 


However, there are many inconsistencies; e.g., 


n several cases, and the heavy rains in the Ohio Valley 
ud New England lie in areas of minimum vertical motion 
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in others. This would indicate that other parameters 
such as available moisture, warm air advection, thickness 
advection, the extent of over-running of a cold stable air 
mass, etc., should be included in any specification scheme. 


Unfortunately the measurement and data processing of 


these parameters is not as easy as for vertical motion. 

This experiment is similar to that done by Woffinden 
[5], and was repeated only because of the difference in 
the circulation and weather regime. However, the general 
results and conclusions are the same. 
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with the understanding that they have not been published 
or accepted for publication elsewhere. 

Two copies of the manuscript should be submitted. All 
copy, including footnotes, references, tables, and captions 
for figures, should be double spaced with margins of at 
least 1 inch on sides, top, and bottom. Some inked cor- 
rections are acceptable, but pages with major changes 
should be retyped. The style of capitalization, abbrevia- 
tion ete., used in the Review is governed by the rules set 
down in the Government Printing Office Style Manual. 

An abstract should be supplied with all but the briefest 
of papers. It should cover, in as concise a manner as 
possible, the main purpose and subject matter of the paper 
and indicate briefly the principal conclusions. One para- 
graph is usually sufficient. 

Tables should be tvped, each on a separate page, with 
u title provided. They should be numbered consecutively 
in arabic numerals. 

In equations conventional svmbols in accordance with 
the American Standards Association Letter Symbols for 
Meteorology should be 
into the manuscript in longhand, dubious-looking symbols 


used. If equations are written 
should be identified with a penciled note. 

References should be listed on a separate sheet numbered 

in the order in which thes occur in the text; or, if there 

than 10, in alphabetical order according to 

The listing should title, 

source (if the source is a magazine, the volume, number, 


are more 


author. include author, and 
month, vear, and complete page numbers should be listed; 
if the source is a book, the publisher, place and date of 
publication, and the page numbers of the reference). If 
a self-contained publication, the 


reference is made to 


author, title, publisher, place of publication, and date 


should be given. Within the text, references should be 
indicated by arabic numbers in brackets to correspond to 
the numbered list. 

Footnotes should be numbered consecutively in arabie 
numerals and indicated in the text by superscripts. Each 
should be typed at the bottom of the page on which the 
footnote reference occurs. 

Illustrations.—A list of captions for the illustrations 
should be typed (double spaced) on a separate sheet, 
Every illustration should have an explanatory caption, 
Number each illustration in the margin or on the back 
To fit into the Review page, 
reduction to 3% by 9 inches 
Map bases 


outside the image area. 
illustrations must 
(column size) or 7: by 9 inches (page size). 
should show only political and continental boundaries 


take a 


and latitude and longitude lines, unless data are to be 
plotted, when station circles will also be needed. Usually 
the less unnecessary detail in the background the better 
will be the result from the standpoint of clear reproduction. 

Line drawings and graphs should also be uncluttered 
with fine background grids unless the graph demands very 
Lines should be even and black, lettering 


both 


close reading. 


legible, svmbols distinctive, and lettering and 
svmbols large enough to maintain their legibility under 
the necessary reduction. If reproductions of line drawings 
are submitted because of the large size of the original, 
they must be extremely clear. Mimeographed copies or 
“photocopy” types of reproduction are not satisfactory. 
Photographs should be sharp and clear, with a glossy 
surface. Bear in mind that marks from paper clips or 
writing across the back will show up in the reproduction. 
Drawings and photographs should be protected with 


cardboard and mailed flat. 
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